BAHAN KURSUS 
DINAMIKA FLUIDA KOMPUTASIONAL 


MODUL II 


ANALISIS DAN PERANCANGAN AEROFOIL 
KOMPONEN BANYAK DALAM ALIRAN POTENSIAL 


Untuk Staf Baru Kelompok Engineering dan Teknologi, 
Calon Spesialis Aerodinamika 


Oleh 


HADI WINARTO B.E.(Syd), M. Eng. Sc, Ph.D (UNSW) 


DEPUTI BIDANG PENGEMBANGAN DINAMIKA FLUIDA 
PUSAT PENGEMBANGAN METODA TEKNOLOGI DAN PRODUKSI 
PT. INDUSTRI PESAWAT TERBANG NUSANTARA 
Jl. Pajajaran 154 Bandung 40174 


BAB I 
PENDAHULUAN 


Konsep dan teori dasar pergerakan fluida pada umumnya telah dibahas secara 
cukup mendalam di modul I buku-buku kursus Dinamika Fluida Komputasional ini. 
Pergerakan fluida viskos dan mampat telah dibahas secukupnya secara kualitatif. 

Dengan menganggap bahwa fluida bersifat tak-mampat, tak-viskos dan tak- 
rotasional telah ditunjukkan bahwa persamaan atur dinamika fluida yang rumit dapat 
disederhanakan menjadi persamaan Laplace yang bersifat linier. Untuk kasus aliran 2- 
dimensional persamaan Laplace dapat diselesaikan dengan bantuan teori variabel 
kompleks. Telah ditunjukkan bahwa solusi dari persamaan Laplace 2-dimensional adalah 
semua fungsi analitik dalam variabel kompleks. Solusi sederhana ini dikenal sebagai 
solusi elementer atau singgularitas. Beberapa diantara singgularitas tersebut mempunyai 
peran yang sangat penting dan dapat diberi arti fisis yang berguna dalam analisis 
aerodinamik. Singgularitas tersebut adalah sumber/serap, doblet dan vorteks. 

Sifat linieritas persamaan Laplace dapat dimanfaatkan untuk memperoleh solusi 
yang lebih rumit sebagai hasil kombinasi linier solusi-solusi elementer tadi. Sebagai 
contoh fungsi potensial kompleks untuk angin seragam, doblet dan vorteks. Komponen 
real fungsi potensial kompleks adalah fungsi potensial sedangkan fungsi arus adalah 
komponen imajinernya. Garis arus adalah kedudukan titik-titik di medan aliran 
disepanjang mana fungsi arus bernilai tetap. Distribusi kecepatan di medan aliran dapat 
dihitung karena kecepatan adalah gradien dari fungsi potensial. Rumus Bernoulli 
memberikan hubungan langsung antara tekanan dengan kecepatan. Ini berarti bahwa bila 
fungsi potensial kompleks suatu aliran itu diketahui maka kita dapat menghitung semua 
kuantitas yang diinginkan. 

Teori pemetaan konformal dengan variabel kompleks dapat dimanfaatkan untuk 
memetakan lingkaran menjadi bentuk aerofoil. Kenyataan ini dapat digunakan untuk 
secara tidak langsung mengevaluasi bentuk potensial kompleks aliran di sekeliling 
aerofoil, seperti bentuk kurva-kurva garis arus, isopotensial, serta distribusi kecepatan 
dan tekanan, dapat dihitung. 

Masalah yang dihadapi kemudian adalah bagaimana caranya mencari bentuk 
rumus transformasi yang diinginkan. Lingkaran dianggap pada bidang komputasi dimana 
informasi tentang garis arus, distribusi tekanan dan lain sebagainya dapat dihitung secara 
analitis karena fungsi potensial kompleksnya diketahui. Aerofoil dianggap berada di 


bidang fisis karena masalah fisis yang ingin dianalisa adalah aliran potensial disekeliling 
bentuk aerofoil tersebut. Garis arus di bidang komputasi dapat secara langsung dipetakan 
menjadi garis arus di bidang fisis dengan menggunakan rumus transformasi. Distribusi 
kecepatan di bidang fisis adalah fungsi dari kecepatan di bidang komputasi dan turunan 
dari fungsi transformasi. Ini berarti bahwa apabila rumus transformasi bentuk lingkaran 
menjadi bentuk aerofoil yang diinginkan itu dapat ditemukan maka semua informasi 
tentang aliran di sekeliling aerofoil tersebut menjadi diketahui. 

Rumus transformasi Joukowski dan Karman-Trefftz adalah 2 rumus transformasi 
yang terkenal dan berbentuk sederhana. Kenyataan tersebut dapat dimanfaatkan untuk 
analisis aerodinamik bentuk-bentuk aerofoil Joukowski dan Karman-Trefftz. Kelemahan 
dari kedua rumus transformasi tersebut adalah kenyataan bahwa aerofoil Joukowski dan 
Karman-Trefftz itu berbentuk tertentu. Kedua rumus tersebut tidak dapat digunakan 
untuk menganalisa aliran di sekeliling aerofoil berbentuk sembarang. Inilah sebabnya 
mengapa kita harus mencari suatu fungsi transformasi yang bersifat lebih umum dan 
mampu memetakan suatu bentuk aerofoil sembarang menjadi bentuk lingkaran atau 
sebaliknya. 

Salah satu metoda pemetaan aerofoil bentuk sembarang menjadi lingkaran adalah 
dengan menggunakan rumus transformasi berbentuk deret Laurent. Dari pertimbangan 
teoritis dapat ditunjukkan bahwa suatu fungsi analitik kompleks selalu dapat ditaksir 
dengan baik oleh deret Laurent, dengan catatan bahwa jumlah suku deret tersebut cukup 
panjang. Masalah yang dihadapi kemudian adalah bagaimana cara mencari koefisien- 
koefisien deret tersebut. Tentu saja koefisien-koefisien tersebut merupakan fungsi dari 
koordinat titik-titik pada permukaan bentuk aerofoil yang ingin dipetakan. Metoda 
hitungan koefisien deret Laurent ini adalah salah satu pokok bahasan dalam buku ini. 
Metoda tersebut dijelaskan secara terperinci sehingga dapat digunakan untuk merancang 
dan membuat sebuah program komputer yang dapat menghitung koefisien-koefisien 
deret tadi. Untuk tinggal landas, pesawat terbang membutuhkan suatu piranti penambah 
gaya angkat seperti "flap", "slat" dan sebagainya. Ini berarti bahwa aerofoil (sayap 2- 
dimensional ) yang perlu dianalisa biasanya terdiri dari banyak komponen. Aerofoil 
tersebut tentu saja dipetakan menjadi banyak lingkaran dan kita harus mampu 
menyelesaikan masalah aliran potensial di sekeliling banyak lingkaran. Untuk kasus satu 
lingkaran, solusi masalah aliran potensial dapat dihitung dengan mudah. Tetapi 
masalahnya menjadi jauh lebih rumit untuk kasus 2 atau lebih lingkaran. Masalah ini 
dapat diselesaikan dengan bantuan teorema lingkaran Milne-Thomson yang pada 


prinsipnya mudah untuk dimengerti. Dalam prakteknya, penerapan teorema lingkaran 
Milne-Thomson untuk menganalisa aliran potensial disekeliling 2 atau lebih lingkaran 
ternyata menjadi cukup rumit dan hitungannya membutuhkan bantuan komputer. Dalam 
buku ini masalah tersebut dibahas secara terperinci. Pemetaan aerofoil komponen tunggal 
menjadi satu lingkaran dapat dilakukan relatif mudah, tetapi untuk kasus aerofoil multi- 
komponen prosedurnya menjadi jauh lebih rumit. Masalah ini juga dibahas sampai tuntas 
dibuku ini. 

Suatu metoda lain untuk menganalisa suatu aerofoil bentuk sembarang adalah 
metoda panel. Konsep singgularitas diskrit dapat diberlaku-umumkan menjadi konsep 
singgularitas diskrit kontinyu. Dalam metoda panel gangguan yang muncul karena 
keberadaan aerofoil dimedan aliran angin seragam dimodelkan dengan menganggap 
bahwa suatu lembaran singgularitas diletakkan pada permukaan aerofoil. Permukaan 
aerofoil ditaksir oleh suatu kurva linier sepotong-sepotong yang berkesinambungan 
membentuk kurva yang mirip dengan bentuk aerofoil. Masing-masing garis linier tersebut 
sesungguhnya adalah perpotongan antara sebuah bidang datar, atau panel, dengan bidang 
dimana aerofoil berada. Inilah sebabnya mengapa metoda ini disebut metoda panel. 

Konsep metoda panel untuk aerofoil multi-komponen tidak berbeda dari metoda 
panel untuk elemen tunggal, jadi metoda ini sangat mudah dan tepat untuk mempelajari 
karakteristik aerofoil multi-komponen. Konsep singgularitas distribusi kontinyu dan 
prosedur metoda panel dibahas secara tuntas dalam buku ini, cukup untuk diterapkan 
dalam membuat program komputer untuk menganalisa bentuk aerofoil sembarang. 

Analisis adalah proses dimana sebuah bentuk aerofoil diberikan dan kita ingin 
menghitung karakteristik aerodinamik bentuk tersebut. Masalah ini dikenal sebagai 
masalah langsung ("direct problem" ) karena metoda panel dapat langsung diterapkan 
untuk menghitung karakteristik aerodinamik aerofoil yang diberikan. 

Dalam masalah perancangan, kita menentukan karakteristik aerodinamik aerofoil 
yang diinginkan dan harus menghitung bentuk aerofoil tersebut. Masalah ini disebut 
masalah balik ("inverse problem" ) karena metoda panel tidak dapat digunakan secara 
langsung untuk mencari bentuk aerofoil. Prosedur masalah balik dapat dijelaskan secara 
ringkas sebagai berikut. Bentuk aerofoil pertama-tama harus ditebak. Aerofoil tebakan 
ini kemudian dianalisa dan hasilnya dibandingkan dengan karakteristik aerodinamik yang 
diinginkan. Perbedaan dari karakteristik aerodinamik aerofoil tebakan dengan yang 
diinginkan kemudian dipakai sebagai data untuk memperbaiki bentuk aerofoil tebakan. 
Prosedur tersebut diulang berkali-kali dalam suatu proses iterasi sampai akhirnya 


konvergen, yaitu karakteristik aerodinamik aerofoil tebakan menjadi sama dengan yang 
diinginkan dan diberikan sebagai masukan. Masalah utama dari prosedur perancangan 
tersebut diatas adalah bagaimana cara memperbaiki bentuk aerofoil tebakan supaya 
prosedur tersebut pasti konvergen dan membutuhkan jumlah iterasi yang sekecil 
mungkin. 

Dalam perancangan dan pembuatan program komputer para pembaca dianggap 
telah menguasai suatu bahasa program komputer, khususnya bahasa FORTRAN. 
Penerapan metoda variabel kompleks dan matoda panel untuk menganalisa karakteristik 
aerodinamik bentuk-bentuk aerofoil, dalam prakteknya melibatkan beberapa masalah 
numerik seperti interpolasi, menghitung solusi persamaan yang rumit, integral numerik, 
menyelesaikan sebuah sistem persamaan aljabar simultan dan lain sebagainya. Masalah- 
masalah numerik tersebut juga dibahas secara rinci di buku ini, supaya para pembaca 
dapat membuat program komputer yang diinginkan cukup dengan membaca bahan-bahan 
yang cukup untuk digunakan membuat program-program komputer yang relatif 
sederhana dan merupakan awal dari program komputer aerodinamika komputasional 
yang dibutuhkan dalam analisis dan perancangan aerofoil, sayap dan bentuk pesawat 


seutuhnya. 


BAB II 
PERSAMAAN-PERSAMAAN ATUR DINAMIKA FLUIDA 


Pada dasarnya semua pergerakan benda diatur oleh hukum-hukum pergerakan New- 
ton , khususnya hukum Newton ke-dua tentang pergerakan benda yang menyatakan 
bahwa gaya berbanding lurus dengan massa dan percepatan (lihat Newton, 1687). Konsep- 
konsep dasar yang berkaitan dengan pergerakan partikel dapat dibaca di buku-buku teks 
fisika (misalnya : Halliday and Resnick, 1964) dan para pembaca dianggap sudah 
mengetahuinya. 

Karena fluida adalah materi kontinyua yang tanpa batasan atau identitas jelas, maka 
hukum-hukum pergerakan fluida menjadi lebih rumit. Pengertian fisis tentang konsep- 
konsep yang diperlukan dalam pergerakan fluida telah diberikan secara gamblang oleh 
Shapiro (1966) dengan pendekatan populer tanpa dibebani oleh matematika yang rumit. 
Pembahasan tentang sifat-sifat fluida dan konsep pergerakan fluida juga telah diberikan 
secara lebih mendalam dan matematis oleh Batchelor (1974). 

Persamaan atur aliran fluida mampat (compressible) tak viskos pertama kali 
diturunkan oleh Euler (1755) dengan menganalisa pergerakan fluida yang berada dalam 
suatu volume kontrol, yaitu volume kecil berbentuk tetap dan terpancang di ruang medan 
aliran. Fluida yang keluar masuk volume kontrol tentu saja silih berganti, berubah dengan 
waktu dan hanya mempunyai identitas waktu berada dalam volume kontrol. Analisis 
pergerakan fluida sebagai partikel dengan identitas tertentu, yaitu yang terdiri dari massa 
yang sama walaupun volume dan bentuknya dapat berubah, telah dilakukan oleh Lagrange 
(1781), tetapi tidak akan dibahas dalam buku ini. 

Euler tidak mengikut sertakan dampak viskositas pada pergerakan fluida dalam 
analisisnya. Perbaikan analisis Euler dengan memasukkan dampak viskositas pertama kali 
dilakukan oleh Navier (1822) dan kemudian disempurnakan Stokes (1845). Untuk 
menghormati Navier dan Stokes maka persamaan atur dinamika fluida yang bersifat umum 
, yaitu mampat dan visko, sekarang disebut persamaan Navier-Stokes. Persamaan atur 
dinamika fluida mampat tetapi tak viskos disebut persamaan Euler. 

Penurunan persamaan atur dinamika fluida telah diberikan oleh banyak penulis, 
misalnya Oswatitsch (1956), Liepmann and Roshko (1956), Shapiro (1958), Schlichting 
(1960), Yuan (1967), Anderson (1982) dan Schreier (1982), jadi tidak perlu diulang lagi 
di sini. Tetapi untuk memperjelas beberapa konsep yang agak sulit dimengerti, di bawah 
ini akan diberikan penjelasan tentang pergerakan solida sederhana yang diharapka dapat 


membantu mengatasi kesulitan dalam membayangkan pergerakan partikel fluida yang jauh 
lebih rumit. 

Marilah kita bayangkan suatu kubus dari zat padat (misalnya kayu atau besi) yang 
berada dalam suatu pipa berpenampang lintang bujur sangkar (lihat gambar 2.1). Bila pada 
kubus tersebut kita terapkan gaya dorong ke arah kanan pada muka AD maka kubus 
bergerak ke kanan. Gerakan ini ditentang atau dilawan oleh gaya gesek yang beraksi pada 
muka AB dan CD ke arah kiri. Bila gaya dorong yang diterapkan lebih besar dari gaya 
gesek, maka sesuai dengan hukum Newton kedua kubus akan bergerak ke kanan dengan 
percepatan atau akselerasi, a, yang diberikan oleh 

Gaya Netto — Massa x Akselerasi 


dimana gaya netto adalah 
Gaya Netto —- Gaya Dorong - Gaya Gesek 

Untuk kasus yang dibahas, pergerakan kubus disebut gerak translasi. Disamping gerak 
translasi ada juga gerak rotasi yang dapat dibayangkan dengan "memaku" kubus ABCD 
pada titik pusatnya, P, dan kemudian menerapakan gaya puntir untuk memutar kubus. 
Dalam gerak translasi seluruh bagian kubus bergerak dengan kecepatan dan akselerasi 
yang sama dengan untuk titik pusat kubus. Sebaliknya dalam gerak rotasi titik pusat kubus 
tidak bergerak, sedangkan titik-titik lain pada kubus bergerak dalam trayektori berbentuk 
lingkaran relatif terhadap titik pusat P. Pada umumnya pergerakan benda adalah gabungan 
atau kombinasi gerak translasi dan gerak rotasi. Karena solida bersifat kokoh, tidak 
berubah bentuk, maka kedua kurva yang mewakili dinding atas dan dinding bawah pipa 
dimana kubus berada (lihat gambar 2.1) harus berbentuk dua lingkaran konsentris, supaya 
kubus ABCD dapat bergerak dalam gerakan kombinasi translasi dan rotasi. Tetapi 
seandainya kubus terbuat dari solida yang sangat elastis dan dapat berubah bentuk dengan 
mudah, maka kurva dinding atas dan bawah pipa tidak perlu berbentuk lingkaran 
konsentris, bahkan pipa tidak perlu sejajar. Kubus elastis tersebut tentu juga dapat 
melakukan kombinasi gerak translasi dan rotasi dengan trayektori sembarang. Untuk kasus 
ini trayektori titik kubus dalam gerak rotasi relatif terhadap titik pusat P tidak harus 
berbentuk lingkaran. 

Pergerakan partikel fluida dapat dibayangkan seperti pergerakan solida yang sangat 
elastis tadi. Dalam membayangkan pergerakan fluida kita tidak membutuhkan "pipa" 
dimana kubus tadi bergerak. Perlu diingat bahwa kubus partikel fluida litu berdempetan 
dengan partikel-partikel fluida lain yang mengerumuninya baik di depan dan di belakang, 
atau bawah dan atas maupun di kedua sisi sampingnya. Semua partikel fluida, bukan hanya 
kubus yang diamati geraknya, bersifat "sangat elastis" yaitu dapat berubah bentuk dengan 


mudah. Untuk pergerakan fluida kurva-kurva dinding pembatas kubus, seperti dijelaskan 
untuk kasus pergerakan solida kokoh, itu sesungguhnya tidak ada. Namun demikian 
konsep tersebut bermanfaat sebagai alat bantu dalam membayangkan pergerakan fluida 
dalam pikiran kita, dan diberi nama garis arus untuk aliran 2-dimensional atau pembuluh 
arus untuk kasus aliran 3-dimensional. 

Bila aliran fluida bersifat tunak ("steady") maka garis atau pembuluh arus akan 
berbentuk tetap, tidak berubah dengan waktu dan pergerakan partikel fluida dapat dengan 
relatif mudah dibayangkan, yaitu mirip dengan kasus pergerakan partikel solida yang 
sangat elastis (lentur) seperti dijelaskan di atas. Sebaliknya bila aliran bersifat turbulen 
(misalnya asap rokok) maka bentuk garis atau pembuluh arus selalu berubah dengan 
waktu, dan dapat saling bertabrakan dan teraduk sehingga identitasnya menjadi tidak jelas 
lagi. Oleh karena itu uraian tentang aliran turbulen tidak dapat diberikan menggunakan 
konsep garis arus dan membutuhkan bantuan konsep-konsep rerata statistik yang sangat 
rumit. Karena buku ini bersifat pengenalan awal maka masalah aliran turbulen tidak akan 
dibahas lebih lanjut. 

Dengan pengertian yang dapat diperoleh dari penjelasan tersebut di atas, para 
pembaca diharapkan akan dapat lebih mudah memahami proses penurunan persamaan 
Navier-Stokes yang diberikan secara terperinci dalam buku-buku teks seperti Yuan 
(1967), Schreier (1982), Anderson (1985) dan lain sebagainya. Dibawah ini hanya hasil- 
hasil akhir dari proses penurunan persamaan atur dinamika fluida yang dibahas tanpa bukti. 
Pada dasarnya persamaan atur dinamika fluida yang dibutuhkan untuk keperluan analisis 
hanya berjumlah dua, yaitu persamaan kontinyuitas dan persamaan momentum. Persamaan 
kontinyuitas adalah model matematis dari hukum Fisika dasar yang dikenal sebagai hukum 
kekekalan massa, sedangkan persamaan momentum dapat diturunkan dari hukum 
kekekalan momentum. Untuk aliran yang melibatkan kecepatan tinggi, perpindahan panas 
dan/atau reaksi kimia ada beberapa persamaan tambahan yang dibutuhkan dan dapat 
diperoleh dari hukum-hukum dasar Termodinamika dan Kinetika Kimia. Contoh dari 
kasus-kasus ini adalah aliran di sekitar pesawat terbang supersonik, aliran gas-gas hasil 
pembakaran dalam ruang bakar, gas buang roket, aliran di sekitar pesawat angkasa ulang 
alik dan lain sebagainya. Dalam buku ini kita hanya membahas kasus-kasus dimana 
persamaan atur dinamika fluida cukup diberikan oleh hukum-hukum konservasi massa, 
momentum dan energi, ditambah dengan beberapa persamaan status yang mungkin bersifat 
semi empiris. 

Bentuk paling umum dari persamaan-persamaan atur dinamika fluida dalam notasi 
vektor dapat diberikan sebagai berikut : 
Persamaan Kontinyuitas : 


Op 2 
—4V.(pV)-0 
2 (pV) 


Persamaan momentum : 


PIL —-VP- 2Vlu(v. 0 V-(20O)JV 4 vafuvsV) (2.2) 
Persamaan Energi : 
Dh DP 
—S—AYH-HD 2.3 
One ar (2.3) 
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Ox Oy Oy Oz 
Patah LB Un 
Oz Ox 3 Ox Oy Oz 
Persamaan Status untuk Gas Sempurna : 
P-—pRT 2.4) 


Persamaan Status untuk Aliran Isentropik : 


Ea konstanta (2.5) 


p" 
Simbol-simbol yang digunakan mempunyai arti sebagai berikut : 
P : densitas atau massa jenis (kg/m” ) 

P : tekanan atau gaya per satuan luas (Pa atau N/m” ) 
h : entalpi (J/kg) 

@ : laju aliran panas (J/detik atau W) 

T : temperatur (" C atau K) 

HO : fungsi disipasi panas (J/detik atau W) 

: viskositas (kg/m/detik) 

: konstanta gas (8.317 J/g mol/K) 

: vektor kecepatan 


“Er 


V-ux4vyiwz 
&5752 adalah satuan vektor di ketiga arah sumbu x,y, dan z yang saling tegak lurus dalam 
koordinat Kartesian. u,v, dan w adalah komponen-komponen kecepatan. Besaran vektor 
kecepatan atau, singkatnya, kecepatan adalah 


V VI 5 (W h 1W3)2 (m/detik) 
Y nisbah panas jenis, yaitu panas jenis proses tekanan tetap dibagi panas jenis proses 


volume tetap atau Cp/Cy. Bernilai 1,4 untuk udara standar. 


V : operator del atau nabla, yaitu vektor GL 4# 13 4 1) 
Vs : gradien dari skalar s 

V. V: divergen dari vektor V (produk skalar V dan V) 
VxV: rot dari vektor V (produk vektor V danV) 


v : viskositas kinematik (- x/p) (m” /detik) 


Persamaan atur di atas berlaku umum untuk fluida materi kontinyua yang memenuhi 
hukum viskositas Newton, atau fluida Newton. Hukum viskositas Newton menyatakan 
bahwa gaya geser berbanding lurus dengan laju kerenggangan fluida dan konstanta 
perbandingannya adalah viskositas. Untuk kebanyakan fluida, khususnya udara, viskositas 
adalah fungsi lemah dari temperatur saja dan tidak tergantung pada tekanan atau sifat 
fluida lainnya. Oleh karena itu untuk aliran fluida yang tidak melibatkan perubahan 
temperatur yang besar, viskositas boleh dianggap konstan dan persamaan (2.2) dapat 
disederhanakan menjadi 


Persamaan Momentum dengan Viskositas konstan : 


DV 2 & 
p—— VP 4 UV'V 4—1V(V.V) 
Dt 3 (2.6) 


Perhatikan bahwa x adalah turunan partikel atau turunan terhadap waktu untuk partikel 


fluida, yang memberikan laju perubahan gerak suatu partikel fluida sebagai fungsi waktu 
dan letak dimana partikel berada, yaitu 


D ORI an 
—-—4V.V 
Dt at (2.7) 


Suku 2 adalah laju perubahan dengan waktu, sedangkan suku V.V adalah laju 


perubahan dengan letak dalam ruang sebagai akibat pergerakan dan dikenal sebagai 
turunan arus atau turunan konvektif. Turunan arus ini perlu diikutsertakan dalam 
persamaan (2.2) karena persamaan tersebut adalah hasil penerapan hukum Newton kedua 


untuk fluida yang berada dalam suatu volume kontrol dan bukan untuk partikel fluida 
dengan identitas atau massa tertentu (massa kontrol). Percepatan gerak bukan hanya 
terbatas pada laju perubahan kecepatan dengan waktu, tetapi juga karena adanya 
perubahan kecepatan dengan letak dalam ruang. Pada waktu tertentu partikel fluida berada 
di sebuah volume kontrol tertentu dan sesaat kemudian berada di volume kontrol yang 
lain. Kecepatan fluida di kedua volume kontrol itu boleh berbeda karena kecepatan adalah 
fungsi ruang disamping fungsi waktu. Ini berarti bahwa partikel fluida mengalami 
perubahan kecepatan karena dia terbawa hanyut dikonveksiokan oleh aliran ke tempat lain. 
Perubahan inilah yang disebut percepatan konvektif atau percepatan arus dan diberikan 
sebagai V.V dalam persamaan (2.7). 

Dengan memanfaatkan definisi (2.7) sekarang persamaan (2.6) dapat diuraikan lebih lanjut 


menjadi 


MA VNV) VP AV VA TA.) 23) 


Sistem persamaan Navier-Stokes yang dijelaskan di atas itu sangat sulit untuk diselesaikan, 
jadi perlu disederhanakan. Untuk aliran dengan kecepatan subsonik rendah, yaitu jauh 
lebih kecil daripada kecepatan suara, fluida dapat dianggap bersifat tak-mampat 
(“incompressible"), yaitu densitasnya berharga tetap. Untuk kasus ini persamaan 
kontinyuitas (2.1) disederhanakan menjadi 

V.V-0 (2.9) 
sedangkan persamaan momentum (2.8) disederhanakan menjadi 











Ng. yy--Iyptovey : (2.10) 
p 
Untuk aliran 2-dimensional, persamaan (2.9) dan (2.10) dapat diuraikan, ditulis dalam 
bentuk lengkap sebagai 
Ea (2.11) 
Ox Oy 
2 z 
AN ea Ha eta pl hi (2.12a) 
ot Ox Oy p Ox Ox Oy 
2 2 
LN men Ka an ata (2.12b) 
Ot” Ox Oy po De Oy 


Persamaan kontinyuitas (2.11) mempunyai banyak solusi dan secara umum dapat 
dikatakan bahwa solusinya adalah suatu fungsi skalar ruang yang disebut fungsi arus dan 
diberi simbol w (x, y). Fungsi arus harus memiliki sifat-sifat berikut 





MEA dan ya (2.13) 
oy Ox 
sehingga 
2 2 
BA 5 eN dan Sa & EA» 
Ox  Ox0y Oy  0yOx 


Karena urutan turunan parsial itu tidak penting, maka diperoleh hubungan berikut 


du WOy yo 
Ox Oy Oxoy OyOx 


yaitu fungsi W memang adalah solusi persamaan (2.11). Suatu kuantitas yang berkaitan 








erat dengan gerak rotasi fluida adalah vortisitas dan vektor vortisitas @, yang dapat 
dirumuskan sebagai berikut : 


O —tx4nytCz (2.14) 
dimana : 
. WN 
Ea 2 (2.15a) 
. du Ow 
biara & (2.15b) 
LD 
Ox Oy (2.15c) 


perhatikan bahwa vektor vortisitas @, sesungguhnya adalah rot dari vektor V, yaitu 
produk vektor dari operator del V, dengan vektor kecepatan, yaitu V, yaitu 

@ srot V-Vx V (2.16) 
Produk vektor dari vektor A dan B diberikan oleh rumus 

A — alif # a27 # a3z 

B -— biz #b27 # b3z 


x y 2 
AxB — lal a2 a3|- (a2b3—a3b2)x #(a3bl—a1b3)y #(alb2—a2b1)z 
b1 b2 b3 
(2.17) 
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Ini berarti bahwa 





x Y 2 
rot V -V x V- 1 “ S3 
Ox Oy Oz 
Uu v w 
Oow Ov 0u Ov Ou. 
“Gay Pa ha ja ra aa YAI (ax aj) (2.18) 


Perbandingan antara (2.14) dan (2.18) menunjukkan bahwa persamaan (2.16) itu memang 
benar. Untuk kasus aliran 2-dimensional di bidang x-y, vektor vortisitas d disederhanakan 








menjadi 
ov Ou 
aa (—x -—)i 
oy 
atau dalam variabel fungsi arus W, sebagai 
C- ov 0u nen, 
Ox Oy Oy 
Yaitu 
Vy —-t (2.19) 


Sekarang persamaan momentum (2.12) dapat diolah lebih lanjut sebagai berikut. 
Persamaan (2.12a) diambil turunan parsialnya terhadap y, sedangkan (2.12b) diambil 


turunan parsialnya terhadap x dengan hasil sebagai berikut 
Ata oU du , PL beat Dv Ou | TE 
ot ay” Ay Ox 2 cp Oy ay. “ay Oy 








2.1 Op | 23 
p Dry” ehi 2 9y: 
d On, Du Ov Dn, WN 0. 


ea Oa OK Ox) OK Oy Dy Dx 





A4 Op, LA LA C (w 
p Bay." Pap 2x) taat 
Selisih dari kedua persamaan di atas setelah diolah lebih lanjut dengan sedikit operasi 
aljabar dan menggunakan persamaan (2.15c) untuk definisi g , akhirnya memberikan hasil 


berikut 
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0G du Ov 0G... Gs 06 0G 
4 4 
at G( dy )»-u—-4-v v—-t—I 








Pemasukan persamaan kontinyuitas (2.11) dan definisi fungsi arus (2.13) ke persamaan di 
atas akhirnya memberikan rumus persamaan momentum dalam variabel fungsi arus, W dan 


vortisitas, c, sebagai berikut : 
Oc 09 06 Op 0g 
Ot OyOx Oxoy 
Seperti dapat dilihat, persamaan Navier-Stokes untuk aliran tak-mampat, 2- 





vVc (2.20) 


dimensional ternyata dapat disederhanakan menjadi persamaan (2.19) dan (2.20) dalam 2 
variabel bebas, w dan g. Persamaan momentum (2.12) yang terdiri dari 2 komponen 
dapat diringkas menjadi 1 persamaan, yaitu (2.20), yang dikenal sebagai persamaan 
transportasi vortisitas. Perlu diperhatikan bahwa persamaan (2.19) dan (2.20) adalah 
persamaan atur untuk fluida Newton tak-mampat dan 2-dimensional. Sistem persamaan 
yang terdiri dari 2 persamaan diferensial parsial, derajat 2 dan tak-linier tersebut masih 
sangat sulit diselesaikan dan perlu disederhanakan lebih lanjut. Persamaan atur dinamika 
fluida itu sulit diselesaikan karena sifatnya yang sangat tak-linier dan ini disebabkan oleh 
suku-suku yang mengandung viskositas. Walaupun dampak viskositas itu sangat penting, 
tetapi ada banyak kasus aliran fluida dimana pengaruh viskositas hanya dirasakan disuatu 
lapisan tipis dekat permukaan benda padat, yaitu lapisan batas. Untuk tujuan 
menyederhanakan persamaan atur dinamika fluida maka dalam analisis yang diberikan 
dibawah, kita akan menganggap bahwa viskositas berharga nol, yaitu dampak viskositas 
akan diabaikan. Apabila viskositas boleh diabaikan maka sistem persamaan Navier-Stokes 
(2.1), (2.2), dan (2.3) dapat disederhanakan menjadi sistem persamaan Euler sebagai 
berikut : 
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Persamaan Kontinyuitas 


2 y.(pV)-0 (2.21) 
Ot 
Persamaan Momentum 
DV DV. a 
——p(—4V.VV)—-—VP 2.22 
Pa ta ) (2.22) 
Persamaan Energi 
Dh Oh - OP - 
—SP(—A4V-Vh) -—4#V-VP 4G 2.23 
Pig Pa ar” B aa 


Walaupun telah disederhanakan, persamaan Euler diatas masih bersifat sangat tak-linier 
dan tidak diketahui solusi analitisnya. Solusi persamaan Euler dapat dihitung secara 
numerik dan merupakan persamaan yang paling sering dipelajari dalam bidang 
aerodinamika komputasional. Analisis dan perancangan aerodinamika pesawat terbang 
transonik biasanya dilakukan dengan menyelesaikan sistem persamaan Euler tersebut. 
Karena buku ini bersifat pengenalan awal maka penyelesaian persamaan Euler tidak akan 
dibahas disini. 

Untuk masalah aerodinamika subsonik, densitas dapat dianggap berharga tetap dan 
persamaan kontinyuitas disederhanakan menjadi bentuk persamaan (2.9). Selanjutnya 
apabila aliran dianggap bersifat tunak maka persamaan momentum dapat diintegralkan 
menjadi Bernoulli sebagai berikut : 


P- hpv? — konstanta (2.24) 


di sepanjang garis arus. 

Untuk lebih memudahkan permasalahan, sekarang dibuat anggapan bahwa aliran 
bersifat 2-dimensional dan tanpa gerak rotasi atau gerak berputar (tanpa pusaran).Aliran 
tanpa pusaran adalah aliran dimana vortisitas berharga nol diseluruh medan aliran. 
Persyaratan atau kondisi tanpa pusaran (irrotasionality condition) dapat ditulis sebagai 
berikut : 


“ Ou 
-. Lg 2.25 
c 2x Oy (2.25) 
Dari persamaan (2.16), persamaan diatas dapat juga ditulis sebagai berikut : 
VxV-0 (2.26) 
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Persamaan (2.26) akan terpenuhi apabila vektor kecepatan V adalah gradien dari suatu 
fungsi skalar ruang, &(x,y), yang disebut fungsi potensial, karena 

V xV6-0 (2.27) 
untuk semua fungsi skalar ruang & 
Ini berarti bahwa untuk aliran tanpa pusaran atau tak-rotasional, vektor kecepatan harus 
memenuhi syarat : 


V — vb (2.28) 
atau 
u-— 2 dan v- 2 (2.29) 
Ox Oy 
Persamaan kontinyuitas (2.11) dapat ditulis dalam variabel & sebagai berikut : 
ek 00 
Pa ay? - V4(x,y)-0 (2.30) 


Persyaratan tanpa pusaran (2.26) juga berlaku untuk kasus aliran 3-dimensional, jadi untuk 
aliran 3-d persamaan kontinyuitas dapat ditulis sebagai berikut : 
ob 00, 0d 
—L-V9(x,y,z2)50 2.31 
Pa ay? ta h(x,y,z) — (2.31) 


Persamaan Laplace (2.30) atau (2.31) adalah persamaan diferensial parsial derajat 2 yang 








bersifat linier dan solusi-solusi elementernya dapat dicari dengan menggunakan pendekatan 
fungsi Green (lihat misalnya Duchateau and Zachmann (1986)) 

Untuk kasus aliran 2-dimensional, solusi juga dapat diperoleh dengan menggunakan 
pendekatan metoda variabel kompleks. Karena buku ini bersifat pengenalan awal maka 
hanya metoda variabel kompleks saja yang akan dibahas secara mendalam. 

Kondisi tanpa pusaran (2.25) juga dapat ditulis dalam variabel fungsi arus, diturunkan dari 
persamaan (2.19), yang disederhanakan menjadi : 


o Ya O'y 
Ox “ap 
Dari persamaan (2.13) dan (2.29) diperoleh hubungan berikut : 








-V'y(xsy)-0 (2.32) 


BG LP gan Pen PM 
Ox Oy oy Ox 
Dalam teori variabel kompleks, hubungan (2.33) dikenal sebagai kondisi Cauchy-Riemann 
(lihat misalnya Churchill (1960), Kreyszig (1962)). Salah satu kesimpulan yang dapat 


(2.33) 


15 


diambil dari hubungan (2.33) adalah bahwa & dan w adalah konyugat kompleks dari 
masing-masing dan suatu fungsi potensial kompleks @(z) dapat dibentuk sebagai berikut : 


O(2)-&(x,y) tiy(x,y) (2.34) 
z—x-tiy 
i-/4 

Ini berarti bahwa persamaan (2.31) dan (2.32) dapat digabung atau disatukan menjadi 
V'9(z)-0 (2.35) 


Persamaan Laplace (2.35) adalah persamaan atur dinamika fluida yang paling sederhana 
dan berlaku untuk aliran tunak, tak-mampat, tak-viskos, tak-rotasional dan 2-dimensional. 
Persamaan tersebut berlaku untuk aliran udara disekitar aerofoil pesawat terbang yang 
melaju dengan kecepatan rendah, kecuali untuk suatu kawasan tipis dekat permukaan 
aerofoil dimana dampak viskositas berperan sangat penting. 

Dalam bab-bab selanjutnya kita akan membahas cara-cara penyelesaian persamaan 
Laplace, baik secara analitis (metoda variabel kompleks), semi-analitis (metoda panel) 
ataupun secara numerik (metoda taksiran beda berhingga). Tetapi sebelum membahas 
masalah penyelesaian persamaan Laplace, ada beberapa konsep dasar dalam bidang 
aerodinamika yang perlu dibahas untuk memberikan pengertian fisis yang lebih mendalam. 
Metoda komputasi walaupun sangat dibutuhkan dan bermanfaat, tetapi sesungguhnya 
hanyalah alat bantu. Dalam masalah aerodinamika yang dibutuhkan adalah pengertian yang 
mendalam mengenai fenomena fisis yang: diamati atau dipelajari. Pengertian tersebut 
kemudian diterapkan untuk menyelesaikan masalah perancangan aerodinamika yang 
dihadapi. 


BAB III 
BEBERAPA KONSEP DASAR AERODINAMIKA 


Aerodinamika pesawat terbang adalah ilmu yang mempelajari bentuk luar pesawat 
terbang dan interaksinya dengan udara yang mengalir di sekelilingnya, dengan tujuan 
menemukan bentuk yang paling tepat untuk memenuhi persyaratan-persyaratan 
perancangan aerodinamika pesawat terbang. Gaya-gaya aerodinamika beraksi pada 
pesawat karena adanya interaksi antara pesawat yang bergerak dengan atmosfir yang 
stagnan (diam), atau antara pesawat yang diam dengan udara yang bergerak seragam, 
seperti dalam terowongan angin. Ada dua.jenis gaya yang beraksi pada pesawat, yang 
memegang peranan kunci dalam perancangan pesawat, yaitu gaya angkat (lift) dan gaya 
tahan (drag). Supaya pesawat dapat dikendalikan naik-turun atau belok kiri-kanan, maka 
pesawat harus dilengkapi dengan bentuk permukaan-permukaan kendali yaitu lap, 
aileron, rudder, elevator dan lain sebagainya pada sayap dan ekor pesawat. Permukaan- 
permukaan kendali tersebut dapat menimbulkan adanya momen-momen aerodinamika 
yaitu momen angguk (Pitching moment), .momen geleng (yawing moment) dan momen 
guling (rolling moment) ' 

Gaya angkat harus dihasilkan oleh bentuk pesawat untuk mengatasi gaya tarik bumi 
(gravitasi) atau berat pesawat yang sedang mengangkasa. Gaya angkat dihasilkan oleh 
sayap dan sirip datar ekor pesawat, dan harus sama besar dengan berat pesawat sehingga 
pesawat berada dalam keadaan kesetimbangan gaya pada arah vertikal (lihat gambar 3.1). 
Gaya angkat ekor datar itu sesungguhnya sangat kecil dan hampir seluruh berat pesawat 
dilawan oleh gaya angkat sayap. Titik dimana gaya angkat sayap beraksi itu tidak tepat 
sama dengan pusat gravitasi pesawat, dan ini berarti bahwa ada momen angguk yang 
cenderung untuk memutar hidung pesawat ke atas atau ke bawah relatif terhadap letak 
ekornya. Oleh karena itu gaya angkat ekor dibutuhkan untuk menghasilkan momen 
angguk yang menyetimbangkan momen angguk tersebut di atas. Dalam buku ini kita tidak 
akan membahas stabilitas atau kendali pesawat dan oleh karena itu gaya yang beraksi pada 
ekor tidak akan dibahas lebih lanjut. Hanya gaya-gaya yang beraksi pada sayap saja yang 
akan dibahas secara terperinci. 

Gaya tahan adalah gaya pada arah yang berlawanan dengan arah pergerakan pesawat 
dan harus diatasi atau disetimbangkan oleh gaya dorong sistim propulsi pesawat. Bentuk 
pesawat harus dirancang sedemikian rupa sehingga gaya angkat yang dihasilkan menjadi 
sebesar mungkin, sedangkan gaya tahan harus diusahakan menjadi sekecil mungkin. Pada 
dasarnya gaya tahan itu muncul karena adanya gelombang kejut dalam aliran mampat 
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(Compressible flow), atau sebagai dampak dari kenyataan bahwa udara sebenarnya bersifat 
viskos. Untuk menyederhanakan masalah, dalam analisis pergerakan fluida yang telah di 
bahas dalam bab sebelumnya kita telah membuat anggapan bahwa fluida yang diselidiki itu 
bersifat tak-mampat, tak-viskous dan tak-rotasional. Anggapan tersebut memang dapat 
menyederhanakan sistem persamaan Navier-Stokes yang rumit menjadi persamaan Laplace 
yang cukup sederhana dan dapat dianalisa lebih lanjut. Tetapi tentu saja hasil analisis tidak 
akan bermanfaat untuk menghitung gaya tahan karena telah mengabaikan dampak 
ketermampatan dan viskositas. P 

Penampang lintang atau profil sayap disebut aerofoil. Pola aliran udara di sekeliling 
sebuah bentuk aerofoil tertentu ditayangkan dalam Gambar.3.1 untuk memberikan 
gambaran tentang apa yang disebut dengan pola aliran. Kurva-kurva di sekeliling aerofoil 
yang digunakan untuk memberikan gambaran tentang pola aliran udara disebut garis-garis 
arus, dan berkaitan erat dengan konsep fungsi arus. Garis arus adalah trajektori atau jejak 
pergerakan sebuah partikel fluida yang mengalir, dalam contoh di atas yaitu dari kiri ke 
kanan. 

Udara dan fluida pada umumnya mempunyai suatu sifat yang disebut tekanan dan 
terdiri dari dua komponen, yaitu tekanan statik dan tekanan dinamik. Tekanan statik 
adalah sifat fluida yang tidak tergantung pada geraknya, sedangkan tekanan dinamik 
berkaitan erat dengan pergerakan fluida. Jumlah tekanan statik dan tekanan dinamik 
disebut tekanan total, yang besarannya konstan atau tetap, tidak berubah di sepanjang 
garis arus. Tekanan dinamik ditentukan oleh kecepatan gerak partikel fluida dan besarnya 
ditentukan oleh rumus berikut 

P, — 2, pv? (3.1) 

2 A 
dimana : 

P, adalah tekanan dinamik (N/m? atau Pa) 

p adalah densitas atau massa jenis (kg/m' ) 
dan V adalah kecepatan (m/detik) 1 
Tekanan statik, tekanan dinamik dan tekanan total fluida dihubungkan oleh rumus berikut: 

P, — P4 P, (3.2) 

di mana P, adalah tekanan total (Pa) dan Pg adalah tekanan statik (P2). 

Karena tekanan total bernilai tetap di sepanjang garis arus maka tekanan statik akan 
mengecil jika tekanan dinamik meningkat dan sebaliknya tekanan statik akan meningkat 
jika tekanan dinamik mengecil. Selanjutnya karena tekanan dinamik berbanding lurus 
dengan kwadrat kecepatan, maka tekanan statik akan mengecil jika kecepatan fluida 


18 


meningkat. Kenyataan ini digunakan untuk menjelaskan proses terjadinya gaya angkat 
sebagai berikut. 

Bentuk aerofoil harus sedemikian rupa sehingga aliran udara yang menyapunya 
bergerak dengan cepat di sepanjang permukaan atas aerofoil dan tekanan statik di situ 
menjadi bernilai kecil. Sebaliknya di sepanjang permukaan bawah aerofoil fluida dibuat 
bergerak lebih lambat sehingga tekanan statik di bagian bawah lebih besar daripada 
tekanan statik di bagian atas aerofoil. Perbedaan tekanan statik antara bagian atas dan 
bagian bawah menimbulkan suatu gaya netto yang beraksi ke arah atas yang disebut gaya 
angkat. Distribusi tekanan di sepanjang permukaan dan gaya angkat yang beraksi pada 
aerofoil jelas sangat berkaitan erat dengan bentuk aerofoil. Inilah inti dari masalah yang 
akan dibahas dalam buku ini. 

Persamaan atur untuk pergerakan fluida tunak, tak-mampat, tak-viskos dan bersifat 2- 
dimensional adalah 

0u du 1 Op 


4 - 
2 ax "ay p Ox 





Man nanah 83) 

Aliran 2-dimensional adalah aliran dalam ruang 3-dimensional di mana sifat-sifat aliran 
udaranya hanya bervariasi di 2 arah sumbu, misalnya x dan y, dan tidak tidak tergantung 
pada arah sumbu ketiga, misalnya z. Pipa pada Gambar.3.2 dapat dianggap sebagai pipa 2- 
dimensional karena fluida yang mengalir dari kiri ke kanan hanya tergantung pada arah x 
dan y, dan sama sekali tidak tergantung pada sumbu z. Pipa bundar pada Gambar.3.3 juga 
bersifat 2-dimensional karena bersifat simetri sumbu (axisymmetric), yaitu sifat aliran 
hanya bervariasi arah sumbu x dan arah radial atau arah sumbu r saja. Apabila jari-jari r 
bernilai kecil, sifat aliran dapat dianggap tidak bervariasi dengan arah radial dan hanya 
tergantung pada arah axial atau arah sumbu x saja. Arah aliran dimana-mana dalam pipa 
dianggap sejajar dengan sumbu x dan oleh karena itu komponen kecepatan v dalam 
persamaan (3.3) berharga nol, sedangkan komponen u berharga tetap di semua titik pada 
penampang lintang pipa, yaitu u hanya bervariasi dengan x. 

Untuk aliran 1-dimensional persamaan (3.3) dapat disederhanakan menjadi 

uSu  1dp (3.4) 


dx p dx 
yang dapat diintegralkan menjadi 
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P4 Lpu 2 - konstanta — P total g (3.5) 


Persamaan (3.5) dikenal sebagai persamaan Bernoulli dan telah diberikan sebelumnya 
sebagai persamaan (3.1) dan (3.2). 

Persamaan Bernoulli tetap berlaku untuk kasus di mana sumbu pipa tidak berbentuk 
garis lurus, tetapi kurva lengkung dengan catatan bahwa sifat-sifat aliran pada penampang 
lintang pipa boleh dianggap tetap atau aliran bersifat I-dimensional (lihat gbr 3.3). Pipa 
kecil atau pembuluh lengkung tersebut lebih dikenal sebagai pembuluh arus dan bila luas 
penampang lintang diperkecil mendekati nol maka pembuluh arus menjadi garis arus. 
Pembuluh arus tidak perlu harus berwujud sebagai pipa nyata, tetapi juga boleh bersifat 
konsep pipa dalam pikiran kita, sebagai alat bantu untuk menjelaskan pergerakan atau 
aliran fluida. 

Laju aliran massa fluida yang melintasi sebuah penampang lintang pipa diberi simbol 


m, yaitu 
m -pVA (3.5) 
di mana : m adalah laju aliran massa . (kg/detik) 
p adalah densitas (kg/m3) 
V adalah kecepatan aliran (m/detik) 
A adalah luas penampang pipa (m2) 


Hukum kekekalan massa menuntut bahwa m berharga tetap di sepanjang garis atau 
pembuluh arus. Karena densitas dianggap tetap, ini berarti bahwa kecepatan aliran V harus 
berbanding terbalik dengan luas penampang A. Pada bagian pembuluh arus yang luas 
penampangnya mengecil, kecepatan aliran harus meningkat. 

Walaupun konsep pembuluh arus telah dijelaskan dengan anggapan aliran bersifat 1- 
dimensional, tetapi sesesungguhnya konsep tersebut berlaku juga untuk aliran 3- 
dimensional. 

Suatu medan aliran 3-dimensional selalu dapat dibayangkan bersifat terbagi menjadi 
sejumlah besar pembuluh arus yang saling menyentuh, berdempet-dempetan. Di bawah ini 
kita akan membahas konsep garis arus untuk aliran 2-dimensional. Marilah kita bayangkan 
suatu aliran udara yang tak-mampat, tunak dan 2-dimensional di bidang x-y. Sebuah titik 
A dengan koordinat (XA,YA) berada di medan aliran dan terpaku di situ, sedangkan titik 
T dengan koordinat (XT,YT) juga berada di medan aliran yang sama tetapi letaknya boleh 
dipindah-pindah (lihat Gambar. 3.4). Kurva ABT dan kurva ACT digambar sembarang 
menghubungkan titik A dan T tersebut. 1 
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Hukum kelanggengan massa menuntut bahwa laju aliran atau fluks massa yang masuk 
melintasi kurva ACT harus sama dengan fluks massa yang keluar melintasi kurva ABT. Ini 
berarti bahwa fluks massa yang melintasi kurva apapun yang menghubungkan A dengan T 
hanya tergantung pada letak T saja karena titik A tidak boleh dipindah. 

Kalau fluks massa diwakili oleh w, maka w hanya tergantung pada atau merupakan 





fungsi skalar dari koordinat (x,y) saja, yaitu 

v5 y(xy) (8.7) 
Fungsi di atas disebut fungsi arus dan hanya mempunyai besaran, tanpa arah, jadi bersifat 
skalar, bukan vektor. Karena w adalah fluks massa maka satuan besarannya adalah 
kg/detik. Kalau fungsi arus di titik A berharga WA dan A dipilih sebagai titik acuan maka 
untuk memudahkan pembahasan, harga WA boleh ditentukan sebagai nol. Fungsi arus di 
titik T diberi simbol WT dan harganya tergantung pada letak T di medan aliran, relatif 
terhadap letak titik A. Dengan perjanjian tersebut, fluks massa antara A dan T adalah wT 
dan harga fungsi arus di sebuah titik adalah ukuran fluks massa yang melintasi kurva 
sembarang yang menghubungkan titik tersebut dengan sebuah titik acuan yang telah 
ditentukan sebelumnya. 

Apabila dalam medan aliran ada 2 titik, yaitu T dan P, disamping titik acuan A (lihat 
Gambar. 3.5), maka fluks massa yang melintasi kurva AP tentu saja harus sama dengan 
jumlah fluks massa yang melewati AT dan TP, yaitu 

YP-yTtyPT : 
atau ditulis ulang menjadi : 

YPT-yT-yP (3.8) 
di mana WPT adalah fluks di P relatif terhadap fluks di titik T. 

Apabila WPT — 0 atau WP - WT maka kurva TP harus berimpit dengan garis arus 
yang melewati titik T, sehingga tidak ada fluks massa yang melintasi kurva TP (lihat 
Gambar. 3.5). Ini berarti bahwa garis arus adalah letak kedudukan titik-titik dalam medan 
aliran di sepanjang mana fungsi arus bernilai tetap. Dengan perkataan lain fungsi arus 
bernilai tetapdi sepanjang garis arus dan fluks massa tidak boleh melintasi garis arus, jadi 
garis arus bersifat seolah-olah mewakili permukaan benda padat yang kedap fluida, Karena 
fluida tidak boleh melintasi garis arus, maka arah vektor kecepatan harus sama dengan 
arah garis singgung pada kurva garis arus. $ 

Kalau As adalah panjang suatu elemen yang sangat kecil pada kurva AT, dengan 
komponen Ax dan Ay pada arah x dan Y, maka fluks yang melintasi As adalah hasil kali 
panjang As dengan komponen vektor kecepatan V yang normal atau tegak lurus pada As, 
yaitu 
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Ay-(V"i)As 


Fluks massa yang masuk melintasi As tersebut kemudian keluar melewati komponen 
panjang Ay dan Ax, masing-masing dengan kecepatan u dan -v (lihat Gambar. 3.5), jadi : 


Ay —uAy-vAx 
Ini berarti bahwa 
Ay 5 (V.a) As —uAy- vAx 


Bila As diperkecil dalam limit mendekati nol, maka kuantitas Aw, As, Ax dan Ay 
menjadi diferensial dw, ds, dx, dan dy dan 


dy -(V.a) ds — udy - vdx (3.9) 


Karena w adalah fungsi skalar dari koordinat (x,y), maka dari teori kalkulus multi-variabel 
kita dapat memperoleh hubungan berikut: 


oy oy 
dy - —dy t — dx 3.10 
Yen Oa (3-10) 
Perbandingan antara (3.9) dan (3.10) menunjukkan bahwa 
Oyw oy 
Pan ra dan --— 3.11 
"—9y Ox Ba) 


Perhatikan bahwa hubungan (3.11) itu tepat sama dengan hubungan (2.13) yang telah 
dibahas sebelumnya di bab 2. Pembahasan di atas telah memberikan arti fisis dari kuantitas 
fungsi arus yang telah diturunkan secara matematis dalam bab sebelumnya. 

Walaupun konsep fungsi arus telah diuraikan dengan membuat anggapan bahwa aliran 
bersifat tunak, tak-mampat, dan 2-dimensional, sesungguhnya konsep fungsi arus bersifat 
lebih umum dan berlaku untuk aliran 3-dimensional, mampat dan bahkan untuk aliran 
turbulen. 

Sekarang marilah kita bahas pergerakan sebuah partikel fluida secara lebih terperinci. 
Pada saat T partikel berada di titik dengan koordinat Ff , yaitu 

F5X9.X ty9.Y t20.2 (3.12) 
Partikel bergerak dengan kecepatan V, dimana 

Vu tv tw. (8.13) 
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Sesaat kemudian, yaitu pada waktu (t # dt), partikel berada di titik (F-- dr ) di mana 
vektor dr adalah 
dr — dx.£ #dy$ tdz2 (3.14) 
dan kecepatan adalah (V # dV ), dimana 
bee iga AN aa (3.15) 
oy Oz 


S|& 


dv (adat Da ay Hap z 


ov ov ov £ 
#(— dx —dyt — d 
ea x ay akn z) y 
He kt Tya Minje 
Ox Oy Oz 


Komponen arah X dari dV dapat ditulis ulang menjadi bentuk berikut 


hn t Da, dy- Da 5 suku I 4 suku II 4 suku III 
Ox Oy Oz 

di mana 
sukuT — SPM 


1|, Ou Oow ov 06u 
kau 55 (La Layar La Ld 
suku Ne Adat (3 -ay | 


Oy 
suka m— Hc kas - Ga )- dz - (D - 2 ay | 


Komponen dV arah y dan z juga dapat diolah dengan cara yang sama. Selanjutnya dV 
dapat ditulis ulang menjadi bentuk berikut 


dV - suku I4 suku II suku III (3.16) 


di mana : 
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suku TM gx-k4 OM. dy Ba 

suku II si “SA d Na 
AL) d Kana pan us 
Ah a 4 j3 )- dx # a 4 Deyay ja 

suku DT — ak Ga - Ga )- dx - (& z Da. ja 
AG Aa par 2. an js 
sil Pe Bpag Ay .ay ja 


Pergerakan partikel fluida pada umumnya dapat dikatakan terdiri dari 2 jenis pergerakan 
dasar, yaitu gerak translasi dan gerak berputar atau rotasi. Di samping itu karena partikel 
fluida bersifat sangat lentur mudah sekali berubah bentuk, maka dalam perjalanannya 
partikel fluida juga pasti mengalami deformasi atau perubahan bentuk. Deformasi yang 
terjadi pada partikel fluida dapat dibagi menjadi dua jenis yaitu deformasi linier dan 
deformasi sudut. Untuk memperjelas konsep tersebut, pergerakan sebuah partikel 2- 
dimensional berbentuk bujursangkar ditayangkan pada gambar. 3.6. 

Kecepatan berputar sebuah partikel diberi simbol Vxgt yaitu suatu vektor dengan 
komponen-komponen sebagai berikut 


Vrot“ Vrx 44 Vey 9 #Vep2) (3.17) 
Kecepatan fluida yang berkaitan dengan gerak rotasi diberikan oleh produk vektor antara 


2. ag 
vektor kecepatan rotasi, Vrot, dan @r yaitu vektor letak titik di mana kecepatan ingin 
dihitung, relatif terhadap titik pusat gerakan berputar, yaitu 


Vrotx dr (Vrx &tVry Pt Vr "2 X(dx Rtdy Y tdz-2) 
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Dari Gambar.3.6 dapat dilihat bahwa komponen arah z dari kecepatan rotasi yaitu V,- 
adalah 


V, 5 (Rear Ranu 


v, H2) (3.19) 


atau 


Ox Oy 
Komponen arah x dan y, yaitu V,, dan V,, dapat diperoleh dari sketsa pergerakan fluida 
2-dimensioanl dibidang y-z dan dibidang z-x dengan hasil berikut 


-M ow Ov 

va 1 3 (3.20) 

v, am) 21) 
aan “Gx 


Persamaan (3.18) sekarang dapat diuraikan menjadi 


Veot X dr - He anjaa- (eBay 
21 oz Ox Ox Oy 
MEN 
211ox Oy oy & 


1 f(ow Ov du Ow 


Perhatikan bahwa di atas itu tepat sama dengan rumus untuk suku III dalam persamaan 
(3.16). In berarti bahwa suku III ruas kanan persamaan (3.16) mewakili gerak rotasi 
partikel fluida. Dari gambar 3.6 dapat juga dilihat bahwa deformasi linear untuk partikel 
fluida 2-dimensional dibidang x-y adalah 


Dlinear — Da ak 4 » dyy 
Ox Oy 


Untuk partikel 3-dimensional, rumus deformasi linear menjadi 


Dlinear Da dx H aa dyy 4 2x dzz (3.23) 
Ox Oy Oz 
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Dlinear — 2 dxx La dyy - 2 dzz (3.23) 
Ox Oy Oz & 


Seperti dapat dilihat rumus di atas itu tepat sama dengan suku I ruas kanan persamaan 





(3.16). Ini berarti bahwa suku I tersebut mewakili deformasi linear. Deformasi sudut 
untuk kasus 2-dimensional dapat diperoleh dari gambar 3.6, yaitu 


Dsudut — (Ea) /dy- ( &as) / | 


Dsudut — TE 4 2) 


Untuk kasus 3-dimensional dapat ditunjukkan bahwa 


Daud —H| (Oa Jaye (Jan 
( 2 (2 2") 5 
#5 —4— dx | —4-— Idz (y 
21 ox” ay Oy Ox 
joe Jana (3.24) 


Rumus (3.24) itu tepat sama dengan suku II ruas kanan persamaan (3.16) jadi suku II 
tersebut mewakili deformasi sudut partikel fluida. Dari pembahasan di atas dapat diambil 
kesimpulan bahwa persamaan (3.16) dapat diuraikan berikut 


atau 





rot 


ad — » 
dV — Dlinear 4 Dsudut- V. xdr (3.25) 


dan pergerakan total sebuah partikel fluida diwakili oleh rumus berikut 
Gerak total — Ve Ve x dri Dlinear 4 Dsudut (3.26) 


Jadi pergerakan partikel fluida melibatkan gerak translasi vw gerak rotasi (Ve x dr), 
deformasi linear (Dlinear) dan deformasi sudut (Dsudut). 

Dalam bab 2 kita telah membahas konsep vortisitas (lihat persamaan (2.14)). 
Perbandingan antara persamaan (2.15) dengan persamaan (3.19), (3.20), dan (3.21) 


menunjukkan bahwa kecepatan gerak putar, Vrot, sebuah partikel fluida itu ternyata sama 
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dengan setengah daripada vortisitas partikel tersebut. Seandainya partikel tersebut tiba- 
tiba dibekukan, maka partikel tersebut akan berputar sebagai benda kokoh dengan 


kecepatan putar Vrat dan vortisitasnya adalah DV. Ini berarti bahwa vortisitas sebuah 
partikel fluida dapat dibayangkan sebagai 2 kali lipat kecepatan putar seandainya partikel 
tersebut bersifat kokoh seperti benda padat. Fluida sangat mudah berubah bentuk dan 
tidak dapat menahan gaya geser, betapapun kecilnya, yang beraksi padanya. Bila ada gaya 
geser yang beraksi pada fluida, maka fluida tersebut akan mengalir atau berubah bentuk 
terus-menerus. Ini jelas sangat berbeda dengan benda padat yang mampu menahan gaya 
geser. Benda padat dapat mengalami perubahan bentuk apabila diberi beban gaya, tetapi 
perubahan bentuk yang terjadi biasanya sangat kecil dan bersifat sementara. Bila gaya yang 
beraksi disingkirkan maka benda padat tadi akan kembali ke bentuk semula. Untuk benda 
padat, konsep kecepatan gerak putar dapat dengan mudah dibayangkan. Sebaliknya, 
konsep kecepatan gerak putar untuk fluida itu agak sulit dibayangkan dan oleh karena itu 
untuk menjabarkan gerak putar fluida dibutuhkan konsep vortisitas yang akan dibahas 
secara lebih rinci di bawah ini. 

Pada gambar (3.6) diperlihatkan sepotong kayu berbentuk segi panjang yang diberi 
lubang di titik S. Melalui lubang tersebut kita tancapkan paku pada meja, sehingga paku 
itu menjadi sumbu putar bila kayu tersebut didorong dan berputar sebagai akibatnya. 
Apabila kecepatan putar kayu diberi simbol @ , maka kecepatan linear titik P pada kayu 
yang berjarak r dari titik sumbu S adalah 

V-0-r (3.27) 
dengan arah tegak lurus pada garis SP sesuai arah putaran. Hasil ini dapat diberikan secara 
matematis dalam notasi vektor sebagai berikut 

V- ox : 
dimana v adalah vektor kecepatan linear, @ adalah vektor kecepatan putar dan r adalah 
vektor letak titik P relatif terhadap S. Vektor @ dirumuskan secara matematis sesuai 
dengan aturan tangan kanan, yaitu bila tangan digenggamkan dan ibu jari dibiarkan lurus 
dan arah putaran adalah dari pangkal ke ujung jari-jari tangan, maka @ berharga positif 
dan arahnya adalah dari pangkal ke ujung ibu jari. Jadi apabila putaran terjadi pada 
halaman buku ini, berlawanan arah dengan putaran jarum jam, maka vektor @ positif 
adalah garis tegak lurus pada halaman buku ke arah mata pembaca. 

Dalam masalah-masalah aerodinamika, arah gerak putar fluida yang sering dibahas 
adalah searah putaran jarum jam. Oleh karena itu dalam perjanjian aerodinamika arah 
putaran jarum jam dianggap sebagai arah @ yang positif, yaitu bertentangan dengan 
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perjanjian matematis. Hal ini tidak perlu menimbulkan kerancuan asalkan saja dijelaskan 
definisi mana yang digunakan pada setiap saat, yaitu perjanjian matematis atau aerodinamis 
Aliran fluida dekat pusat atau mata tornado (lihat gambar 3.8) bersifat sebagai gerak 
putar benda kokoh untuk fluida, yaitu distribusi kecepatan fluida dari titik pusat sampai 
jarak tertentu, misalnya r, , diberikan oleh persamaan (3.27). Kecepatan tangensial fluida 
di titik pusat tornado bernilai nol dan mencapai harga maksimalnya, yaitu Vmax, di titik 
berjarak r, dari titik pusat. Kecepatan tersebut kemudian menurun dengan semakin 
meningkatnya jarak radial r dan kecepatan untuk rr, diberikan oleh rumus 


VW) - Vmax. (r, /r) untuk Pr, (3.28) 
Jadi untuk fluida di luar kawasan inti tornado, yaitu untuk rr, , kecepatan tangensial 
berbanding terbalik dengan jarak radial. Sebaliknya di dalam kawasan inti kecepatan 
tangensial berbanding lurus dengan jarak radial, yaitu diberikan oleh rumus berikut 


V(r) - Vmax. (r/r,)  untukrsr, » (3.29) 


Apabila dirumuskan bahwa @, adalah kecepatan gerak putar benda kokoh untuk fluida 
dalam kawasan inti tornado, yaitu 

@,- Vmax/r, (3.30) 
maka persamaan (3.29) dapat ditulis ulang menjadi 

Vr) -0,"r untuk rsSr, (8.31) 
Gerakan-gerakan berputar dalam aliran fluida disebut pusaran. Supaya dapat dibahas lebih 
lanjut dan dipelajari secara kuantitatif, pusaran dibayangkan atau diidealisasikan bersifat 
seperti tornado. Pusaran mempunyai kawasan inti dimana fluida bergerak putar seperti 
benda kokoh sedangkan di luar kawasan inti distribusi kecepatan tangensial berbanding 
terbalik dengan jarak radial dari titik pusat pusaran. Vorteks adalah pusaran dimana jari- 
jari kawasan inti pusaran, yaitu r,, diambil limitnya mendekati nol tetapi harga 
(Vmax“r, ) adalah suatu konstanta yang tidak mendekati nol. Kecepatan tangensial di 
titik pusat vorteks, yaitu Vmax, dianggap tidak terdefinisikan, sedangkan di titik-titik 
lainnya kecepatan tangensial diberikan oleh rumus (3.28) yang ditulis ulang menjadi 


V(r) - K/r (3.32) 


dimana K adalah suatu konstanta, yaitu 
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lim Vmax "r, -K 


Konstanta K disebut kekuatan vorteks dengan satuan besaran m?/detik dan berkaitan erat 
dengan kuantitas lain yang disebut sirkulasi vorteks. Sirkulasi vorteks biasanya diberi 
simbol T dan diberi nilai berikut 


T-21K (3.34) 


Sirkulasi adalah suatu ukuran gerakan berputar aliran fluida pada umumnya dan tidak 
terbatas pada vorteks saja. Konsep sirkulasi yang berlaku umum ini dapat dijelaskan 
sebagai berikut. Dalam suatu aliran 2-dimensional kita dapat menggambar sebuah kurva 
tertutup sembarang dan kita sebut sirkit C. Kurva tersebut adalah suatu alat bantu 
matematis dan tidak mewakili suatu permukaan fisis, jadi sesungguhnya keberadaan sirkit 
C itu hanya ada dalam pikiran kita saja. Vektor kecepatan di titik T pada sirkit C (lihat 
gambar 3.9) adalah V, dan dc adalah panjang sebuah elemen sirkit di T. Sirkulasi 
dirumuskan sebagai integral atau hasil penjumlahan dari komponen kecepatan V yang 
sejajar dengan garis singgung pada sirkit C dikalikan dengan panjang elemen dc, yang 
diambil limitnya mendekati nol, di setiap titik yang berada di sepanjang sirkit C. Dalam 
bahasa matematis sirkulasi I' « pada sirkit C diberikan oleh rumus berikut 

Te-pv-ae-p V cos 9 dc (3.35) 
dimana 0 adalah sudut yang diapit oleh vektor V dan garis singgung pada sirkit C di titik 
T (lihat gambar 3.9). Sirkit C boleh dipilih berbentuk sembarang, jadi untuk lebih 
memperjelas masalah dan memudahkan pembahasan kita pilih suatu bentuk segi empat 
(lihat gambar 3.9). Sirkulasi di sekeliling bentuk segi empat ABCD dapat dihitung sebagai 
berikut 
Rumus untuk vektor V dan dt adalah 

V-u-x4v-y 

de - dx-x4dy-y 
jadi produk skalar kedua vektor tersebut adalah 

V-dt —udx-vdy (3.36) 
dan persamaan (3.35) dapat ditulis ulang menjadi 


Tc- htuax vdy) 
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Penerapan rumus di atas untuk segi panjang ABCD akhirnya memberikan harga sirkulasi 
sebagai berikut 


Txseo “ Uxp-AB 4 V,o.BC-— U.,,. CD — V,,.DA (3.37) 


Kuantitas Uk,» Ucps Vpc dan V,, adalah harga-harga rerata dari kecepatan u dan v 
di sepanjang sisi-sisi AB, CD, BC dan DA dari segi empat ABCD. Apabila ukuran segi 
empat ABCD diperkecil dan diambil limitnya mendekati nol, yaitu AB —- CD — 6x dan BC 
-DA - 6y, dan 6x,8y — 0 maka sirkulasi ABCD menjadi 





£ Vsc— V Um-U 
ST Anc-jim (0x-00| un DA - 2 
0y—0 


Selanjutnya ditentukan bahwa koordinat titik pusat segi empat ABCD adalah (x,y) 
sehingga 


vtxt201)- a10) u(y4205)-u(y—0y) 


ST ABcp-lim (0x-0y) - & 


2y-0 





atau 


Oo oOu 
dr — (dx- dy —— — 
ABCD — (dx: dy) ay) 
Karena elemen segi panjang ABCD berada dibidang x-y dan arah kebalikan putaran jarum 
jam diberikan sebagai satuan vektor arah sumbu-z positif, maka 


df, — dT,-2 (dx- na (3.38) 
Vortisitas di suatu titik T dalam medan aliran dirumuskan sebagai suatu vektor arah tegak 
lurus pada bidang (baik bidang datar atau kurva ruang) dimana suatu sirkit C yang 
berukuran kecil tak-berhingga dapat dibayangkan mengurung titik T tersebut. Harga 
besaran vortisitas adalah limit harga sirkulasi pada sirkit C dibagi oleh luas bidang yang 
terkurung dalam sirkit C apabila luas tersebut diambil limitnya mendekati nol. Karena 


30 


sirkit C diambil limitnya dengan luas mendekati nol maka sirkit tersebut pasti berada pada 
suatu bidang datar 2-dimensional dan arah vektor vortisitas terdefinisikan dengan jelas. 
Vektor vortisitas diberi simbol & karena kemiripannya dengan konsep kecepatan gerak 
putar untuk benda kokoh dan dalam ruang 3-dimensional mempunyai komponen- 
komponen sebagai berikut 
O —txtnytz (3.39) 

Karena vortisitas adalah sirkulasi per satuan luas, maka komponen @ arah z adalah (lihat 
persamaan 3.38) 


—dT, 0 Ou 


dx.dy Ox Oy 


Persamaan (3.38) memang telah diperoleh untuk suatu bentuk sirkit C tertentu, yaitu 
segi empat ABCD, tetapi di bawah ini akan dibuktikan bahwa hasil di atas berlaku umum 
untuk semua bentuk sirkit C. Apapun juga bentuk sirkit C, sirkit ini selalu dapat dibagi 
menjadi sejumlah besar elemen segi empat yang amat kecil, dengan sisi-sisi sejajar dengan 
sumbu-x dan sumbu-y (lihat gambar 3.10). Sirkulasi di kawasan yang dikurung oleh sirkit 
C adalah jumlah dari semua sirkulasi pada elemen-elemen segi empat karena kontribusi 
atau sumbangan dari suatu sisi pembatas dua elemen, yang saling berdempetan, terhadap 
sirkulasi total itu berlawanan tanda dan saling meniadakan, jadi yang tersisa hanyalah 
sirkulasi kontribusi dari kurva pembatas atau sirkit C saja. Ini berarti bahwa untuk setiap 
sirkit C berbentuk apapun, rumus untuk sirkulasinya adalah 


re: Par, 


atau dengan memanfaatkan rumus (3.38) maka 


TA 5 2 Ja dy (3.40) 
Tetapi sirkulasi juga diberikan oleh rumus (3.36), jadi 
@v Ou 
udx - vdy) — ——— Idx.d 3.41 
ih y) H -- 2| y (841) 


dimana A adalah luas yang dikurung oleh sirkit C. Untuk aliran 3-dimensional, sirkulasi 
diberikan oleh rumus berikut 


Ie - pV.dz- fP(o-n)aA | (8.42) 


A 
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dimana 


$V- ae $(udx 4 vay 4 wdz) (8.43) 
c 


c 


fa aan (We Jane M- Manan 4( 2 Janay| 


A A 
(3.44) 
Ini berarti bahwa komponen-komponen vektor vortisitas, @ , adalah 
Ow Ov 
Sa 3.45a 
$ ay” Gx (3.45a) 
Ou Ow 
Se 3.45b 
PA Ag (8.45b) 
-X Ou (3.456) 
Ox Oy . 


Karena vortisitas atau sirkulasi adalah ukuran untuk aliran dengan pusaran atau gerak 
berputar, maka aliran tanpa pusaran adalah aliran dimana vortisitas atau sirkulasi berharga 
nol di setiap titik dalam medan aliran. Aliran tanpa pusaran disebut aliran potensial dan di 
bawah ini konsep aliran potensial akan dibahas secara rinci. Dalam suatu medan aliran 
terdapat dua titik, yaitu A dan T, dimana A dipaku di satu tempat sedangkan T boleh 
dipindah-pindahkan (lihat gambar 3.11). Sirkulasi pada sirkit ABTCA adalah 


pv-ae- | V-act JV.ae (3.46) 
ABCTA ABT TCA 
Mengingat bahwa 
D 5 L 
(v.de—— Jv.dc (3.47) 
TCA ACT : 


maka persamaan (3.46) dapat ditulis menjadi 


pv-ae- | V-ae— |V.ac (3.48) 
ABCTA ABT ACT 


Apabila aliran yang diselidiki adalah aliran potensial maka sirkulasi pada sirkit ABTCA 
pasti berharga nol, jadi 
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(Vde- (V.de- JV-de 69 
AC A 


ABT 


Seperti dapat dilihat integral dari kuantitas (V-dt) untuk aliran potensial itu tidak 
tergantung pada jejak integrasi dan hanya tergantung pada titik-titik ujung limit integral. 
Karena titik A terpaku di satu tempat maka harga integral hanya tergantung pada letak T 
saja. Ini berarti bahwa kita dapat merumuskan suatu fungsi skalar ruang yang harganya 
hanya tergantung pada letak titik dimana harga fungsi dihitung. Fungsi skalar ruang 
tersebut diberi nama fungsi potensial dengan simbol & dengan sifat sebagai berikut 





dp -(V-dc) (3.50) 
sehingga 
T F4 E. 
db, — Jap -— | V-ac (3.51) 
A A 
Mengingat bahwa 
V-de-—u.dx # v.dy #w.dz (3.52) 
sedangkan dari teori kalkulus multi-variabel kita tahu bahwa 
de dea Ciaya Mae (3.53) 
Ox Oy Oz 
maka dapat diambil kesimpulan berikut 
V - Vo (3.54) 
atau 
a30. 
Ox 
u& 
Oy 
“8 
Oz 


Hasil di atas menunjukkan bahwa untuk aliran potensial, vektor kecepatan adalah 
gradien dari fungsi potensial. Kurva yang merupakan letak kedudukan titik-titik dengan 
harga fungsi arus yang sama disebut kurva isopotensial. Karena kecepatan adalah gradien 
dari fungsi potensial maka arah pergerakan fluida adalah normal atau tegak lurus pada 
kurva isopotensial. Tetapi dari konsep garis arus diketahui bahwa arah aliran itu sejajar 
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dengan garis singgung pada garis arus, jadi dapat diambil kesimpulan bahwa kurva 
isopotensial dan garis arus itu saling tegak lurus satu dengan lainnya. Ini sebabnya 
mengapa fungsi potensial 6 dan fungsi arus w dapat dianggap sebagai konyugat 
kompleks, yang selalu saling tegak lurus, dalam teori variabel kompleks. 

Dalam bab ini kita membahas arti fisis dari beberapa konsep dasar dalam bidang 
aerodinamika seperti fungsi arus, fungsi potensial, vortisitas, sirkulasi, persamaan 
Bernoulli dan lain sebagainya. Konsep-konsep tersebut akan dipelajari lebih lanjut dalam 
bab-bab berikut, khususnya dalam kaitan dengan pembahasan mengenai aliran potensial, 
tak-mampat. 


34 


BAB IV 
FUNGSI ANALITIK VARIABEL KOMPLEKS 


Salah satu cabang Matematika adalah ilmu berhitung atau Aritmatika yang 
mempelajari masalah bilangan. Dalam aritmatika pada awalnya kita belajar tentang 
bilangan bulat positif dan operasi matematis yang paling sederhana yaitu penjumlahan, 
atau tentang operator #. Selanjutnya kita juga belajar tentang operasi-operasi matematis 
yang lain, yaitu pengurangan (-), perkalian (x), pembagian (/), perpangkatan, logaritma 
dan lain sebagainya. 

Apabila bilangan yang kita kenal hanya terbatas pada bilangan bulat positif, ternyata 
satu-satunya operasi matematis yang dapat dilakukan adalah penjumlahan. Supaya 
operasi pengurangan dapat dilakukan tanpa batas, kita harus merumuskan bilangan nol 
dan bilangan bulat negatif. Dengan sistem bilangan yang mencakup bilangan bulat positif, 
bilangan nol dan bilangan bulat negatif kita dapat melakukan operasi tambah, kurang dan 
perkalian. Tetapi supaya kita dapat melakukan operasi pembagian, kita harus 
merumuskan bilangan pecahan rasional, misalnya 1/2 Sistem bilangan bulat dan pecahan 
rasional merupakan suatu sistem yang lengkap dipandang dari segi operasi-operasi 
aritmatika dasar yaitu 4, -, x dan /, dan disebut sistem bilangan rasional. Kemudian dalam 
aljabar kita belajar tentang konsep variabel dan persamaan aljabar. Suatu contoh 
sederhana dari persamaan aljabar adalah persamaan kuadrat berikut 

x—450 
dengan solusi-solusinya yaitu x — t2 
Tetapi untuk kasus yang sedikit berbeda, yaitu 

x—2-0 
ternyata sistem bilangan rasional itu tidak cukup lengkap. Untuk menangani masalah di 
atas kita harus merumuskan suatu jenis bilangan baru yaitu bilangan pecahan tak-rasional. 
Salah satu contohnya adalah bilangan V2 yang merupakan solusi dari persamaan 
x — 2-0. Sistem bilangan rasional dan tak-rasional itu ternyata juga belum lengkap dan 
harus ditambah dengan bilangan transendental, misalnya bilangan alamiah, e, dan nisbah 
keliling dibagi diameter lingkaran, t, yang merupakan solusi dari persamaan 
transendental. Sistem bilangan yang mencakup semua jenis bilangan tersebut di atas 
dikenal sebagai sistem bilangan nyata, atau sistem bilangan real, yang dapat diwakili 
secara grafis oleh suatu garis bilangan, yaitu garis bilangan real. Semua titik yang berada 
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pada garis bilangan real selalu dapat diwakili oleh salah satu bilangan real, baik bilangan 
bulat, nol, pecahan rasional atau tak-rasional dan transendental. 

Walaupun sistem bilangan real ini kelihatannya cukup lengkap tetapi ternyata tidak 
mencakup solusi dari suatu persamaan aljabar sederhana berikut 

x”#1-0 
Untuk keperluan ini kita harus merumuskan suatu jenis bilangan baru yaitu bilangan 
khayal atau bilangan imajiner dengan simbol i dan dirumuskan sebagai berikut 

?--1 
atau 

i-J4 
Dengan tambahan bilangan imajiner ini, sistem bilangan kita menjadi lengkap. Semua 
persamaan aljabar sekarang mempunyai solusi atau akar berupa bilangan. Sebagai 
contoh, solusi dari persamaan 

x t4x47-0 
adalah 

x-—24ti/3 
Perhatikan bahwa harga x diatas tidak dapat diwakili oleh suatu titik pada garis bilangan 
real. Dalam kenyataannya, sistem bilangan yang paling lengkap tidak dapat diwakili oleh 
sebuah garis apapun, tetapi harus diwakili oleh titik-titik yang berada di sebuah bidang 2 
dimensional yang dikenal sebagai bidang kompleks. Sistem bilangan yang mencakup 
bilangan real dan bilangan imajiner disebut sistem bilangan kompleks dan diwakili oleh 
titik-titik pada bidang kompleks, yang mempunyai sumbu real dan sumbu imajiner. 
Dengan mengamati kenyataan bahwa 

i-—lxi 

—-ixi 
dan sudut yang dibentuk oleh garis bilangan dari 0 ke 1 dan dari 0 ke -1 adalah 1800 
dapat diambil kesimpulan bahwa perkalian dengan i berarti pemutaran garis bilangan real 
melalui sudut 900. Ini berarti bahwa sumbu bilangan imajiner harus tegak lurus pada 
sumbu bilangan real (lihat gambar (4.1)). Apabila sumbu bilangan real positif ditentukan 
searah dengan sumbu-x positif maka sumbu bilangan imajiner harus sejajar dengan 
sumbu-y dan bidang x-y adalah bidang bilangan kompleks, yang juga disebut bidang-z. 
Bilangan seperti -2 4-i1/3 adalah bilangan kompleks yang diwakili oleh titik pada bidang 
kompleks dengan koordinat (-2,1/3). Pada umumnya suatu variabel kompleks dengan 
simbol z dapat ditulis sebagai x--iy atau sebagai (x,y). Teori-teori tentang variabel dan 
fungsi kompleks telah dibahas oleh banyak penulis (misalnya lihat Churchill (1960), 
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Phillips (1963), Kreyszig (1962) dan lain sebagainya) dan tidak perlu diulang di sini. 
Namun demikian dalam bab ini kita akan membahas teori variabel kompleks secara 
singkat dan secukupnya dalam kaitan dengan cara-cara penyelesaian persamaan Laplace, 
yang dibutuhkan untuk menyelesaikan masalah aliran potensial, yang menjadi inti bahasan 
buku ini. Bilangan kompleks dapat dirumuskan secara formal sebagai berikut. Bilangan 
kompleks adalah suatu bentuk teratur pasangan bilangan real (x,y) yang mematuhi 
hukum-hukum dasar aljabar, yaitu hukum-hukum komutatif dan asosiatif untuk 
penjumlahan dan perkalian, dan juga hukum distributif. Simbol yang digunakan untuk 
bilangan kompleks adalah z, yaitu z-(x,y), sedangkan operasi-operasi aljabarnya dapat 
dijelaskan sebagai berikut 

Bila z, —(x,,y,) dan 2», 5(X5,y2)) » maka 


Z1 2) 5 (Xi ) # Ora ya) Si try, Hy) 

Z7 (Xu Yi DX (X7, Ya) S (X3 — YiYa X4 Ya Hy) 

Zr “2g 5 (X1 Yi) — (X2 y2) 5 Ox yi —y2) (4.1) 
Lia (aya) (aa tYiya XYi | 

Z2. (Xyz) x tyo 5 XP tyo 





Pada prakteknya bilangan kompleks sering ditulis sebagai 
z5(x,y)—xtiy (4.2) 
i-J4 
Komponen x dari bilangan kompleks z disebut komponen real, dengan simbol Re(z), 
sedangkan y disebut komponen imajiner, dengan simbol Im(z). Konyugat kompleks dari 
z-(x,y) adalah bilangan kompleks dengan simbol Z, dimana 7 - ( x,—y), yaitu komponen 
realz dan Z itu sama, sedangkan komponen imajinernya berlawanan tanda, yaitu 


Re(z)-Re(z) | (4.3) 


Im(z) —-Im(z) 
Dalam pembahasan di atas kita telah menggunakan sistem koordinat Kartesian untuk 
menentukan letak titik z di bidang kompleks. Letak titik z dapat juga ditentukan dengan 
menggunkan sistem koordinat polar dan cara ini disebut pemanfaatan diagram Argand 
(lihat gambar (4.2)). Besaran bilangan kompleks diwakili oleh panjang garis dari titik asal 
(0,0) ke titik (x,y) dan disebut modulus bilangan kompleks z, dengan simbol mod(z) atau 
ll, yaitu 


mod(z) — al 5-27 4 y? (4.4) 


37 


Arah garis yang menghubungkan titik asal (0,0) ke titik z disebut argumen, dengan 
simbol arg(z), yaitu 

arg(z) —tan"(y/x) (4.5) 
Dalam diagram Argand atau sistem koordinat polar, bilangan kompleks z ditulis sebagai 
berikut 

z-—(r,0)-r.exp(id) —r(cos0,sin 0) (4.6) 
mana r adalah modulus z dan 0 adalah argumen z. Konyugat kompleks dari z-—(r,0) 
adalah z —(r,-0). Kuantitas r dan @ dihubungkan dengan x dan y oleh rumus berikut 


y-rsin@ (1) 


x5r 3) 
Operasi perkalian dan pembagian dengan bilangan kompleks dapat dilakukan dengan 
lebih mudah menggunakan sistem koordinat polar dibandingkan dengan koordinat 
Kartesian. Sebagai contoh, rumus (4.1b) dan (4.1d) dapat ditulis ulang dalam notasi 
diagram Argand sebagai berikut 

ZXZ, 1,.exp (10, Jr, .exp(i0,) —rr,.exp(i(O, -0,)J 

2 MB)  erpfi(O, 65) 
z, ,.exp(id,) r, 


atau 


zyxz, — (1,,0,)x(r,0,) — (1r,,0, tO,) 

n-tsn.( R00) C3) 

z2 (n,0,) r, 
Apabila x dan y (atau r dan 0) adalah variabel-variabel real, bukan bilangan real, maka z 
adalah variabel kompleks. Kita dapat merumuskan fungsi variabel kompleks, seperti 
halnya dengan untuk variabel real. Fungsi variabel kompleks diberi simbol f(z) dimana z 
2 (x,y) dan dapat diuraikan menjadi 

f2) -f(xy) -g(xy) tih(x,y) — (g,h) (4.9) 
dimana g dan h masing-masing adalah fungsi variabel real x dan y, dan disebut 
komponen real dan komponen imajiner dari fungsi kompleks f, yaitu 


g(x,y)-Relf (z)) | 
h(x,y)- Im(f (z)) 

Supaya lebih jelas kita berikan contoh, yaitu 
f(2)-27 

yang dapat diuraikan menjadi 


(4.10) 
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F2) 5 (xtiyP 5 —y) ti(2xy) 
yaitu untuk kasus ini 

2(x,y) - Retf (z)) —x? —y? 

h(x,y)- Im(f (z)) — 2xy 
Domain dari fungsi f(z) adalah seluruh bidang kompleks z, kecuali titik-titik dimana 2) 
tidak terdefinisikan. Titik dimana harga f(z) tidak terdefinisikan disebut titik singguler 
atau Singgularitas dari f(z). Sebagai contoh, domain dari f(z) —2? adalah seluruh bidang- 
z, sedangkan domain dari f(z) -1/z adalah seluruh bidang-z kecuali titik asal z - (0,0) 
dimana (1/z) tidak terdefinisikan. Titik (0,0) adalah titik singguler atau singgularitas dari 
fungsi (1/z). Fungsi variabel kompleks dapat diturunkan atau diambil integralnya, sama 
seperti untuk variabel real, dengan catatan bahwa harga-harga limitnya itu ada atau 


terdefinisikan. Derivatif atau turunan dari suatu fungsi real, f(x), dirumuskan sebagai 
berikut 





"9-3 im, ME 

Domain bilangan real adalah garis bilangan real dari —co ke --oo, jadi harga limit di atas 
didekati dengan memperkecil harga Ax mendekati nol, menelusuri sepanjang garis 
bilangan tersebut. Tetapi domain bilangan kompleks adalah bidang 2 dimensional, jadi 
harga limit harus didekati dari segala arah, bukan sepanjang satu garis saja, pada bidang 
tersebut. Harga limit suatu fungsi kompleks di suatu titik dapat dihitung dengan 
menggambar suatu lingkaran kecil dengan titik tadi sebagai pusatnya, dan kemudian 
jejari lingkaran diperkecil dalam limit mendekati nol. Turunan dari fungsi kompleks f(z) 
kemudian dapat dirumuskan sebagai berikut 

f (z . f(z24Az)-f(z 

d a Ju. Tim, ( 2 (z) 

dengan catatan bahwa harga limit di atas itu ada. Supaya lebih jelas, di bawah ini kita 





f'(2)— 


berikan satu contoh. Bila f (z) -2?, maka f'(z)— 2z untuk semua z, karena 


lim 
Azb0 


1 De: 
SAN TA EN NA TA 
Az Az1b0 


yaitu f'(z) terdefinisikan untuk semua z. Tidak semua fungsi kompleks f(z) itu 
mempunyai turunan. Sebagai contoh marilah kita periksa fungsi f (2) -Jel, atau 


f(z)-x' 4y?, yang bersifat kontinyu atau mulus dalam variabel real x dan y. 


2 2 
f'(2)- lim ana 
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Kr Ar) 4 (y #AYYI- (247) 








- lim 
Ax,hy-0 Ax #iAy 
yaitu 
#2) Tim (2x Ax)- Ax4-(2y 4 Ay)- Ay 


Ax,Ay-0 -Ax ti- Ay 
Untuk kasus khusus dimana z - (0,0) maka 


2 2 
fta)— Tia KAV AAN) ima (Ar 2iAy)-0 
Ax4y»0 — Ax-#iAy Ax,Ay—0 


Jadi turunan untuk f (z) — lzP itu pada umumnya tidak terdefinisikan dan hanya ada di z — 
(0,0) dengan nilai f"(z)-0. Apabila turunan dari fungsi f(z) itu terdefinisikan, maka 
dapat dipastikan bahwa fungsi tersebut bersifat mulus atau kontinyu. Tetapi sebaliknya 
bila fungsi tersebut bersifat kontinyu, belum tentu bahwa turunan fungsi itu pasti ada, 
seperti telah ditunjukkan dalam contoh-contoh yang diberikan di atas. Supaya turunan 
fungsi f(z) itu pasti ada, maka fungsi tersebut harus memenuhi persyaratan-persyaratan 
tertentu yang akan diselidiki di bawah ini. Pertama-tama kita jelaskan dulu apa yang 
dimaksud dengan pernyataan bahwa f(z) bersifat kontinyu. Suatu fungsi f disebut 
kontinyu di titik z, hanya kalau fungsi tersebut memenuhi tiga persyaratan di bawah ini. 


1. f(z)) terdefinisikan atau ada 
2. limf(z) terdefinisikan atau ada (4.11) 
Z—29 


3. limf(z2)-f(2,) 


Kalau f -g--ih bersifat kontinyu, maka dapat dipastikan bahwa g dan h adalah fungsi 
kontinyu dalam variabel real x dan y. Sebaliknya bila g dan h bersifat kontinyu dalam 
variabel real, maka f pasti juga bersifat kontinyu dalam variabel kompleks. Selanjutnya 
bila fungsi w(z) bersifat kontinyu, maka fungsi komposit f(w(z)), disingkat f(w), pasti 
juga kontinyu. Ini berarti bahwa fungsi dari suatu fungsi kontinyu pasti juga bersifat 
kontinyu. Sebagai contoh, karena g-x? —y? dan h-—2xy adalah fungsi-fungsi kontinyu 
dalam variabel real x dan y, maka fungsi w(z), dimana w(z)-z? -g-4-ih, pasti juga 
kontinyu. Selanjutnya fungsi komposit f(w), misalnya f (w) —-2w tiw? —2z? 4iz', pasti 


juga kontinyu. Sekarang ingat bahwa f (2) le (Pt) rio adalah suatu fungsi 
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kontinyu, tetapi turunannya tidak terdefinisikan. Jadi ada atau tidak adanya turunan suatu 
fungsi ternyata membutuhkan persyaratan yang lebih ketat daripada persyaratan untuk 
sifat kontinyu saja. Kita tentukan bahwa fungsi f (z) — g(x,y) tih(x,y) memiliki turunan 
di Z9 X9 tiyo, yaitu f'(z,)—a—ib, dimana a— a(xo»yo) dan b —b(x,,yp). Sekarang 
marilah kita periksa hubungan antara g, h, x, y, a dan b. Bila ditentukan bahwa 


Af 5f (xp #Az) —£ (2) — Ag-#iAh 
Ag — g(x9 # Ax, yg # Ay) — 8(xX9»yo) (4.12) 
Ah -h(x9 # Ax,y9 #Ay) —h(x9s yo) 
maka harga limit dapat dihitung sebagai berikut 
. Af. AgtiAh 
Tim —— lim ————- 
4150 Az A250 Ax -4- iAy 
Dari teori tentang limit dan pengertian mengenai komponen-komponen real dan imajiner, 


a-ib (4-13) 


kita peroleh persyaratan berikut 


lim Ref Sean) -a 


Ax,Ay—0 Ax #iAy (4.14) 
Tim ia Sean) -—b 
Ax,Ay—»0 Ax #iAy 


Untuk menyederhanakan masalah, kita boleh mengambil kasus khusus yaitu untuk 
kondisi Ay —0, sehingga Az-—Ax dan kedua limit di atas disederhanakan menjadi limit 


fungsi dari satu variabel saja, yaitu 








lim &(Xo#Ax, ya) —8(xosya) -a atau 3 di (X95yYo) 

Ax—0 Ax Ox (4 15) 
“ ilh(xo t Ax, yo) -h(xosyo)I.. Oh 9 
jim, Aa —-ib atau ai b di (x,»yo) 


Proses limit dapat juga diambil sepanjang sumbu-y, yaitu dengan menentukan bahwa 
Ax — 0 sehingga Az—iAy dan 
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ilh(Xo»Yo #Ay) —h(xo»yo)I Oh 











jim, iAy a atau ay -a di (X95Yo) tega 
lim B(Xos Yo t Ay) —8(Xo» Yo) ih tata Sp di (Xosyo) 
Ay-0 iAy Oy 


Jadi persyaratan yang harus dipenuhi, supaya f -g--ih mempunyai turunan di titik 
Z9 5 X9 tiyo, Yaitu f'(z,) —a—ib, adalah adanya hubungan-hubungan berikut 


0g Oh 
Ox 

2 & (417) 
oy 


- 


Sh 
Ox 


Persyaratananya adanya turunan di atas, yaitu hubungan (4.17), dikenal sebagai kondisi 
Cauchy-Riemann dan menjadi dasar teori fungsi analitik dan penyelesaian persamaan 
Laplace. Hasil di atas dapat diringkas menjadi teori yang menyatakan bahwa apabila 
turunan fungsi f(z) itu ada, maka kondisi Cauchy-Riemann pasti dipenuhi, tetapi ini tidak 
berarti bahwa bila kondisi Cauchy-Riemann dipenuhi maka turunan f(z) itu pasti ada. 
Teori kondisi Cauchy-Riemann tersebut dapat ditulis lengkap sebagai berikut Kalau 
turunan dari fungsi f -g ih itu ada di titik z, maka turunan-turunan parsial derajat 
pertama relatif terhadap x dan y dari setiap komponen f, yaitu g dan h, itu pasti ada di 
titik z tersebut dan memenuhi kondisi Cauchy-Riemann, persamaan (4.17). Persyaratan 
yang lebih ketat untuk keberadaan turunan f(z) diberikan oleh teori keberadaan turunan, 
yaitu Bila g dan h adalah fungsi-fungsi real nilai tunggal (single valued) dari variabel 
real x dan y, dan fungsi-fungsi tersebut beserta turunan-turunan parsialnya adalah 
Jungsi kontinyu di titik (x,,yo) , dan kalau turunan parsial tersebut memenuhi kondisi 
Cauchy-Riemann di titik itu, maka dapat dipastikan bahwa turunan f'(z,) dari fungsi 
f —g-ih itu ada, dengan persyaratan bahwa z-x 4iy dan 1, -x, tiy,- Dalam sistem 
koordinat polar, kondisi Cauchy-Riemann dapat ditulis sebagai berikut 


2g Adit 


Or rod 
108. Oh 


r 09 Or 


(4.18) 
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dimana : z-(r,0 ), f(z) —g(r,@) tih(r,0) dan r # 0 

Teori keberadaan turunan dalam sistem koordinat polar dinyatakan sebagai berikut. 
Ditentukan bahwa g(r,@) dan h(r,@ ) masing-masing mempunyai satu harga real disetiap 
titik z di kawasan tetangga sekitar titik (r,,0,). Ditentukan juga bahwa g, h beserta, 
turunan-turunan parsial derajat pertama dari g dan h relatif terhadap r atau @ ,adalah 
fungsi-fungsi kontinyu dari variabel z di (r,,0,). Bila kedua persyaratan diatas dan bila 
kondisi Cauchy-Riemamn dititik (r,,8,) dipenuhi, maka dapat dipastikan bahwa turunan 
f'(z) itu ada, dengan catatan bahwa f -g-#ih, z,-r,.exp(i@,) dan z-r.exp(id) 
sedangkan f'(z) diberikan oleh rumus 


FG) -(cos0 —isin0( B4 MJ -ept-i0) SJ toy 
aja teen iatanl 0 mr il 0 
#'(2)- (cos0 isinO)( 2 128 )--expt io) ig) 


Apabila turunan f'(z) itu ada, bukan hanya di z, tetapi juga disetiap titik dalam kawasan 
tetangga di titik z,, maka fungsi f(z) disebut fungsi analitik. Fungsi f(z) disebut analitik 
disuatu domain pada bidang-z, kalau fungsi itu analitik disetiap titik dalam domain 
tersebut. Selanjutnya suatu fungsi disebut analitik total bila fungsi itu bersifat analitik di 
setiap titik pada bidang-z. Bila fungsi f(z) bersifat analitik untuk semua titik di kawasan 
tetangga titik z,, kecuali di titik z,, maka z, disebut titik singguler atau z, adalah 
singgularitas fungsi tersebut. Kawasan tetangga titik z, adalah kawasan yang dibatasi 
oleh suatu lingkaran kecil, berjari-jari bukan nol, dengan titik z, sebagai pusatnya. 

Bila f(z) dan w(z) adalah fungsi-fungsi analitik, maka fungsi jumlah (f-w), selisih (f- 
w) dan hasil kali (fx w) pasti juga bersifat analitik di domain yang sama, sedangkan hasil 
bagi (f/w) juga analitik di domain yang sama kecuali di titik-titik dimana w(z) berharga 
nol.Selanjutnya bila w(z) adalah suatu fungsi analitik maka fungsi komposit f(w) pasti 
juga analitik di domain yang sama. 

Sekarang marilah kita pilih suatu fungsi f(z) khusus, yaitu fungsi potensial kompleks 
dengan simbol &(z). Komponen-komponen real dan imajiner dari & diberi simbol 
6(x,y) dan w(x,y), yaitu 

D(2z) -&(x,y) tiy(x,y) (4.20) 
Apabila fungsi O(z) bersifat analitik, maka turunannya pasti ada dan komponen- 
komponen & dan w harus memenuhi kondisi Cauchy-Riemann berikut 
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Da 02 


Ox” Oy 
9. 
Oy Ox 


Turunan-turunan parsial derajat kedua dari 6 dan w adalah sebagai berikut 











Oo 09. Oy . 0. 

Ox? Oyox Ox 

dp. 0y. v 

Oy”— Ox0y Oy 
dan juga 

Op Ow. Ou 

Ox0y  Oy”  Oy 


CK MN LIA MUKA. 


OyOx Ox? Ox 
Sekarang perhatikan bahwa 








du, v Ba Nu Oy Owy. 
Ia Oh Dai ay” OyOx OxOy 


dan 





dd V Owy ,0v ad op Ow 
Oy Ox Oy” Ox  Oxdy OyOx 
Persamaan (4.22) dan (4.23) dapat ditulis ulang sebagai 
Vp-0 
V'y-0 
atau digabung menjadi 
V?0(2) 5 V (diy) - Wp tiV'y 0 





-0 


(4.21) 


(4.22) 


(4.23) 


(4.24) 
(4.25) 


(4.26) 


Persamaan (4.26) menunjukkan bahwa apabila O(z) bersifat analitik maka O(z) adalah 


solusi persamaan Laplace. Dalam sistem koordinat polar operator Laplasian, yaitu V?, 


diberikan oleh rumus berikut 
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? 10 1? 
LNG aa 4.27 
OP rar og nd 
dan persamaan Laplace untuk & dan @ adalah 


v 





Lah 1 06 

r Or t r? gg 1 

2y loy 10w Pena 
4——-4 5 

Ore ror ri 99 


V0) St 





V'y(r,0) - 





Kondisi Cauchy-Riemann dalam sistem koordinat polar (lihat persamaan (4.18)) dapat 


ditulis sebagai berikut 


Ob .10y .y 
LL (4.29) 
Loh ...8W 
r 00 Or : 
sehingga 


op. 1 oy 1 Oy.OV, 

AA ee ar un 
Or r' 00 roddr Or (4.30) 
1906 19y 10, 


2 





rdrod rr? 900” r@ 


lOh 196 OwW.ONg 


PE t- 
rd rodor Or Or (431) 
196. 19Pw  10Vp 


r? Oo? roro r & 
Sekarang perhatikan bahwa 








V, 10Vo 06,106 10y 190y 10 
Or ro Or re 99 rr 0 rai rirad 


Persamaan di atas dapat diolah lebih lanjut menjadi 














2 2 
29,196, LOW, 
Ol ro ra 
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dan karena 


ep Ie 

r' 0 rOor 
maka akhirnya 

24 lap 1 F4 

LET SEPENAA ANN AA PET 

Or roro? 
Selanjutnya perhatikan juga bahwa 





1V, VW 196 aa Op 106 10y,0wy 
rod Or rirdd r' 00 rodar re 90? Ar 
Persamaan di atas dapat diolah lebih lanjut menjadi 








Ow, low 1 Ob 


or ra ra 





dan karena 
A0h. A0w 
r? 09 r Or 


maka akhirnya 


2 2 
0 ha slow 4# 10w —0 
Or or or ro? 
Seperti dapat dilihat, apabila 46 dan w memenuhi kondisi Cauchy-Riemann dalam 





koordinat polar (persamaan (4.29)) maka dan w adalah solusi dari persamaan Laplace 


dalam koordinat polar (persamaan (4.28)). Turunan pertama dari fungsi potensial 
kompleks, &(z), dalam sistem koordinat Kartesian adalah 


D' (z)-u-iv (4.32) 
dan konyugatnya adalah 

D' (2) -u-tiv (4.33) 
Dalam sistem koordinat polar rumus untuk turunan dari fungsi potensial kompleks adalah 

D'(2)5(V, —iV,)(cos0 —isin0) (4.34) 
dan konyugatnya adalah 

D' (2) 5(V, tiV,)(cos0 tisin0) (4.35) 


Bila @(z)-&-iy adalah suatu fungsi analitik kompleks maka komponen- 
komponennya, yaitu 4 dan w, disebut fungsi-fungsi konyugat harmonik. Fungsi & 
adalah konyugat harmonik dari fungsi w dan sebaliknya w adalah konyugat harmonik 
dari 6 .Bila salah satu dari kedua fungsi konyugat harmonik tersebut diketahui rumus 
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analitiknya, maka fungsi yang lain dapat dicari rumusnya dengan menerapkan kondisi 
Cauchy-Riemann. Sebagai contoh, bila rumus untuk & adalah 














9-y -3xy 
dan fungsi ini memenuhi persamaan Laplace karena 
Ob ob 2 2 
—5-6 dan ——3y' — 3x 
2x xy Zy y 
ob op 
ana --by dan ay - by 
sehingga 
@ 
v9 21,28 .9 
Ox Oy 


maka rumus untuk fungsi konyugat harmoik dari &, yaitu fungsi w, dapat diperoleh 
dengan mudah sebagai berikut: 
Kondisi Chauchi-Riemann adalah: 


NY dc 

ay Ex 

Dg 23 ag 
Ox Oy 


Integral dari kedua persamaan tersebut adalah 


y- ja 2 Jay - z -| Gxydy — —3xy” 4 F(x) konstanta 


w- ( V Jas sit 3x? —3y 2)ax—x? —3xy? 4 Gly) konstan ta 


Ini berarti bahwa 
—3xy' #F(x) 5x —3xy'4Gly) 
jadi 
Fx) —x dan G(y)-0 
Oleh karena itu rumus untuk wadalah 
W-x' — 3xy' #konstanta 
Dalam bab ini kita telah menunjukkan bahwa solusi dari persamaan Laplace 2- 


dimensional adalah semua fungsi analitik kompleks. Persyaratan yang harus dipenuhi oleh 
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fungsi analitik kompleks juga telah dibahas secara rinci. Beberapa contoh dari fungsi 
analitik elementer (sederhana) adalah 


AG Ta | 1 : 
2 pasal Ia asin(z),exp(z),sin(z),.cosh(z) 


z 
dan lain sebagainya. Mengingat bahwa persamaan Laplace itu bersifat linear, maka kita 
dapat memperoleh solusi-solusi yang lebih rumit, yang memenuhi kondisi-kondisi batas 
tertentu, dengan menggabungkan dalam suatu kombinasi linear beberapa (2 atau lebih) 
solusi-solusi elementer. Persamaan Laplace dan solusi-solusinya telah dibahas dari segi 
pendekatan murni dan belum dikaitkan dengan masalah aliran potensial. Namun demikian 
dari bab-bab sebelumnya kita tahu bahwa model matematis dari aliran potensial tak- 
mampat adalah persamaan Laplace. Oleh karena itu bahan-bahan yang telah dibahas di 
sini akan dimanfaatkan untuk bab selanjutnya untuk mencari solusi masalah aliran 
potensial. 
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BAB V 
POLA-POLA ALIRAN POTENSIAL ELEMENTER 


Dari bab-bab sebelumnya diketahui bahwa persamaan atur pergerakan fluida tunak, 
tak-mampat, tak-viskos, tak-rotasional dan 2-dimensional adalah : 


Persamaan Kontinyuitas 
Ou 3 KA —0 
Ox Oy 
Kondisi ketak-rotasionalan (kondisi tanpa pusaran) 
KM 
Ox Oy 


Disamping itu persamaan momentum dapat diintegralkan dan memberikan persamaan 
Bernoulli berikut 


P- Lpv? — konstan ta 


yang berlaku di sepanjang garis arus dan 

V4 
Supaya hasil-hasil analisis berlaku umum, maka semua kuantitas seperti u, v, p, x, y dan 
lain sebagainya perlu dibuat menjadi kuantitas tanpa dimensi. Ini dapat dilakukan dengan 
menentukan parameter-parameter acuan, seperti kecepatan acuan, panjang acuan, tekanan 
acuan dan seterusnya. Untuk aliran di sekeliling bentuk aerofoil harga-harga acuan adalah 
harga parameter dalam aliran bebas yang tidak terganggu oleh keberadaan aerofoil, 
misalnya kecepatan dan tekanan arus bebas. Acuan untuk panjang adalah panjang aerofoil. 
Sekarang kita dapat merumuskan parameter-parameter tanpa dimensi sebagai berikut : 
Kecepatan tanpa dimensi 


sa... 
u -—— 
Va 
aa 
v -—— 
Va 
vN 
V, 


& 


Tekanan tanpa dimensi 
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P' - P 
Lv 2 
2 Ya 





Panjang tanpa dimensi 


Kuantitas V,, dan c adalah kecepatan dan panjang acuan, sedangkan TpV," adalah 


tekanan acuan. Informasi tentang tekanan biasanya diberikan dalam bentuk koefisien 
tekanan yang dirumuskan sebagai berikut : 











Ca 1 (EP 51- (VP 
TA Vo 


Persamaan kontinyuitas dan kondisi tanpa pusaran dalam variabel tanpa dimensi sekarang 
dapat ditulis sebagai berikut : 














Persamaan Kontinyuitas : 
Ba. bg 2 2 
Ox Oy 

Kondisi tanpa pusaran 
Ov Ou - 
3 


Seperti dapat dilihat, bentuk kedua persamaan dalam variabel berdimensi dan tanpa 
dimensi itu ternyata sama. Untuk memudahkan penulisan simbol bintang, “, tidak akan 
ditulis, tetapi tetap dengan pengertian bahwa semua variabel yang terlibat adalah variabel 
tanpa dimensi. Seperti telah dijelaskan dalam bab-bab sebelumnya, persamaan kontinyuitas 
dan kondisi tanpa pusaran itu dapat digabung menjadi satu, yaitu menjadi persamaan 
Laplace dalam variabel fungsi potensial kompleks, D(z), yaitu 
V9) 5 V9 HiV2y —0 
Komponen real, &, disebut fungsi potensial kecepatan, sedangkan komponen imajiner, w, 


disebut fungsi arus dan hubungannya dengan kecepatan u dan v adalah sebagai berikut 
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na te 


aa 

Ns SE 2 
Oy Ox 

V jo ai 


Semua fungsi analitik kompleks adalah solusi dari persamaan Laplace dan dalam bab ini 
kita akan membahas 4 fungsi analitik elementer secara rinci. 
Keempat fungsi elementer tersebut adalah : 


1. &(z) —z.exp(-—ia) aliran angin seragam (5.1) 

2. O@(z)- toln(z) aliran sumber/serap (5.2) 

3. @(z) — tiyIn(z) aliran vorteks (5.3) 

4. O(z) — B aliran doblet (5.4) 
z 


Pola-pola aliran yang diwakili oleh keempat solusi elementer tersebut diatas dan arti 
fisisnya akan dibahas lebih lanjut dibawah ini. 

Karena @(z)-6-tiw dan z—-x-#iy, maka @(z)—z.exp(—ia) dapat diuraikan 
sebagai berikut 


b(x,y) —xcos(a) -ysin(o) 


y(x,y) -ycos(a) — xsin(a) (5.5) 
Komponen kecepatan u dan v untuk kasus ini adalah 
u(x,y) — cos(a) 3 


v(x,y)- sin(a) 
Karena kecepatan aliran mempunyai besaran dan arah yang tidak berubah dengan 

koordinat (x,y), arti fisis dari solusi ini adalah aliran angin seragam dengan kecepatan V-1 
dan arahnya membentuk sudut « dengan sumbu-x. Trayektori atau jejak pergerakan 
partikel fluida disebut garis arus, yang secara matematis dapat dirumuskan sebagai letak 
kedudukan titik-titik di medan solusi yang mempunyai nilai fungsi arus yang sama. Untuk 
kasus ini rumus matematis untuk garis arus adalah sebagai berikut 

ycos(a) — xsin(a) — w, (kostan ta) 
atau diolah lebih lanjut menjadi 

ysxtan(a) -w, / cos(a) (5.7) 
dimana &« dan w, adalah konstanta. 
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Jelas bahwa pola aliran untuk kasus ini adalah garis-garis arus yang berbentuk garis- 
garis lurus yang membentuk sudut « dengan sumbu-x (lihat gambar (5.1)). Kurva-kurva 
isopotensial diberikan oleh rumus berikut 





Pe AA (5.8) 
tanf(a) sin(a) 


jadi kurva-kurva isopotensial dan garis arus saling berpotongan tegak lurus satu dengan 
lainnya. Solusi elementer kedua yang diwakili oleh persamaan (5.2) dapat diuraikan 
menjadi 


6(x,y) In J2 4y? (5.9) 
w(x,y) —0tan” (y/x) 
Komponen kecepatan untuk kasus ini adalah 





x 
u(x,y) Or 1 
i (5-10) 
y 
xy) -—0——— 
v(x,y) Tg 
Komponen kecepatan dalam koordinat polar adalah 
o $ 
V, (7,90) — — (arah radial) 
r 
Vo (r,0)—0 (arah tangensial) (5.11) 
dimana r — /x? -y? dan@ —tan”(y/x) 
Rumus matematis untuk garis arus adalah 
otan” (y/x)—-00 -— w, (konstanta) (5.12) 
sedangkan kurva isopotensial diberikan oleh rumus 
oln/x #y?' —0Inr — 6, (konstanta) (5.13) 
Rumus (5.12) dan (5.13) dapat ditulis ulang sebagai 
0-w,/o(konstanta) (5.14) 
r — exp(b, /0)(konstanta) (5.15) 


Untuk kasus ini garis-garis arus adalah garis-garis radial yang memancar dari atau ke 
titik asal sistem koordinat. Bila & berharga posistif maka aliran memancar keluar dari titik 
asal dan kasus ini disebut aliran sumber. Sebaliknya bila d berharga negatif maka aliran 
memancar masuk atau menuju ke titik asal dan disebut aliran serap. Kecepatan aliran 
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berbanding terbalik dengan jarak radial. Pola aliran untuk kasus ini dapat dilihat pada 
gambar (5.2). 

Aliran sumber secara fisis dapat dibayangkan sebagai berikut. Bayangkan 2 plat datar 
berbentuk bujur sangkar dan berukuran sangat besar, mendekati tak-berhingga. Kedua plat 
tersebut diletakkan berdekatan dan sejajar satu diatas yang lain. Rongga atau jarak 
pemisah kedua plat tersebut dibuat cukup kecil. Pada plat bawah dibuat sebuah lubang 
bundar kecil dan disambungkan kesuatu pipa yang berada dibawah plat (lihat gambar 
(5.3)). Pipa ini dihubungkan ke suatu pompa air yang bekerja memaksa air dalam pipa 
bergerak keatas dan akhirnya memancar keluar memasuki rongga antara plat. Air akan 
memancar ke rongga, ke segala arah dengan kekuatan pancar yang sama, tidak bergantung 
arah. Karena aliran ini mirip suatu sumber air, maka disebut aliran sumber. Seandainya 
pada awalnya rongga plat itu sudah berisi air dan arah pemompaan dibalik, maka air dalam 
rongga akan tersedot atau terserap dari segala arah masuk kedalam pipa dan aliran ini 
disebut aliran serap. Parameter & mempunyai arti yang sama dengan w, yaitu fluks massa 
tanpa dimensi dan tentu saja menentukan kekuatan pancar aliran sumber. Oleh karena itu 
Oo disebut kekuatan sumber, bernilai positif untuk aliran sumber dan negatif untuk aliran 
serap. 

Suatu lingkaran berjari-jari r dengan titik pusat berada di titik sumber dapat 
dibayangkan digambar pada plat datar tadi. Fluks massa yang memancar dari titik sumber 
adalah w dan sesuai hukum kekekalan massa nilainya harus sama dengan fluks massa yang 
melintasi lingkaran r tadi. Bila kecepatan radial dari fluida yang melintasi lingkaran diberi 
simbol Vy, maka nilai fluks massa adalah 

Y-— 21r9V, 
dimana & adalah rongga atau jarak pemisah kedua plat datar. Ini berarti bahwa kecepatan 
radial diberikan oleh rumus berikut 

1 

v, - j3 (5.16) 
yaitu V, berbanding terbalik dengan r, seperti diberikan oleh rumus (5.11). 
Apabila harga r diambil limitnya mendekati nol, maka kecepatan radial V, menjadi tak- 
terdefinisikan, yaitu r — 0 adalah titik singgular. Walaupun dalam pembahasan diatas kita 
hanya membayangkan 2 plat datar, sesungguhnya kita dapat saja membayangkan sejumlah 
besar plat yang saling sejajar dan terpisahkan oleh rongga setinggi 9 satu dari yang lain. Ini 
berarti bahwa sesungguhnya aliran ini adalah aliran 3-dimensional. Namun demikian pola 
aliran di semua rongga, yang sejajar dengan bidang x-y, itu tepat sama dan tidak 
tergantung pada arah sumbu-z, jadi dianggap bersifat 2-dimensional. Oleh karena itu aliran 
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ini dapat diwakili oleh suatu garis sumber atau garis serap yang sejajar dengan sumbu-z 
dan panjangnya tak-berhingga. Perpotongan garis sumber dengan bidang x-y adalah 
sebuah titik, tetapi disebut garis sumber karena titik tersebut sesungguhnya adalah 
proyeksi garis sumber pada bidang x-y. Hal ini perlu ditegaskan disini karena disamping 
garis sumber ada juga konsep titik sumber yang bersifat 3-dimensional, yaitu fluida 
memancar dari satu titik ke semua arah diruang 3-dimensional (lihat Gbr (5.4)). Karena 
aliran ini bersifat 2-dimensional, harga & dalam persamaan (5.16) boleh diganti dengan 1 
satuan, jadi 


.y1 


Na 21r 


Perbandingan dengan rumus (5.11) menunjukkan bahwa apabila d mempunyai arti fisis 
sebagai kekuatan sumber atau fluks massa maka seharusnya & — w dan rumus (5.2) untuk 


aliran sumber seharusnya diberikan oleh 


O(2) - —inz (5.17) 
2x 

Demikian juga aliran vorteks dan doblet seharusnya diberikan oleh rumus-rumus berikut 
ANE 

P(z) - —i—Inz (5.18) 
21 
H1 

02) Te (5.19) 
21z 


Aliran elementer ketiga yang akan dibahas adalah aliran vorteks yang fungsi 
potensial kompleksnya diberikan oleh persamaan (5.18). Fungsi potensial kecepatan dan 
fungsi arus untuk kasus ini diberikan oleh rumus-rumus berikut 


b(x,y) - tan (S 2 


w(sy) - “ag V4 (5.20) 
T 


Komponen kecepatan untuk aliran vorteks adalah 





Na San 
uG5y)- sa 4y 
x 
vG,y) - ET ag (5.21) 


Komponen kecepatan dalam koordinat polar adalah 
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V. (r,0)-0 








1 (5.22) 
Vetr,0)— -—— 
21r 
Rumus matematis untuk garis arus adalah 
Yk 
. 3 5.23 
r er 2, A (5.23) 
sedangkan rumus untuk kurva isopotensial adalah 
Pe (5.24) 
y/2x 


dimana y, w, dan 6p adalah konstanta 

Ini berarti bahwa kurva garis arus untuk aliran vorteks berbentuk lingkaran-lingkaran 
dengan jari-jari diberikan oleh persamaan (5.23) sedangkan kurva isopotensial adalah 
garis-garis radial yang memancarkan dari titik singgularitas vorteks. Pola aliran vorteks 
dapat dilihat pada gambar (5.5). Kurva isopotensial aliran vorteks ternyata sama dengan 
garis arus untuk aliran sumber, dan sebaliknya garis arus aliran vorteks ternyata sama 
dengan kurva isopotensial untuk aliran sumber. Oleh karena itu dapat dikatakan bahwa 
aliran vorteks adalah konyugat harmonik dari aliran sumber. Suatu contoh aliran vorteks 
yang sering muncul secara alamiah adalah tornado atau angin puting beliung. Suatu contoh 
lain yang dapat dengan mudah diciptakan di laboratorium adalah aliran air dalam suatu 
bejana yang sangat besar dan berlubang kecil pada dasar bejana. Pada awalnya lubang 
tersebut disumbat dan bejana diisi air sampai penuh dan dibiarkan sampai air dalam bejana 
menjadi tenang, diam tak bergerak. Setelah itu sumbat pada lubang di dasar bejana 
kemudian dicabut dan air dibiarkan bocor mengalir keluar bejana. Tidak lama kemudian 
kita akan melihat gerakan air berputar dengan diameter yang cukup besar pada permukaan 
air. Pusaran air ini turun ke bawah menuju lubang di dasar bejana dan diameter pusaran 
tersebut makin mengecil semakin mendekati lubang. Bila bejana terbuat dari plastik atau 
kaca jernih yang tertembus cahaya, maka kita dapat melihat melalui sisi bejana suatu aliran 
vorteks yang mirip dengan tornado, tetapi dalam skala kecil (lihat gambar (5.6)). 

Rumus untuk fungsi potensial kompleks (5.18) telah diberikan dalam perjanjian 
matematis, sehingga apabila y bernilai positif maka arah kecepatan tangensial adalah 
kebalikan arah jarum jam (lihat persamaan (5.22) dan gambar (5.5)). Dalam perjanjian 
aerodinamika, rumus untuk fungsi potensial kompleks aliran vorteks adalah 


@(z) -i--in(z) (5.25) 
2x 
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sehingga untuk harga y positif arah kecepatan tangensial adalah searah putaran jarum jam. 
Fluks massa yang melintasi garis radial dari jari-jari r, sampai ke jari-jari r, dapat 
dihitung sebagai berikut 


. ?ydr y r, 
- 'vdr— (-—-4n2 5.26 
Wrur, | 6ar 15 r 21 n( ) ( ) 


Seperti dapat dilihat, y disini tidak sama dengan w yaitu y bukanlah ukuran fluks massa 
untuk aliran vorteks. Arti fisis dari kekuatan vorteks dapat dipelajari lebih lanjut dengan 
menghitung sirkulasi di sekeliling lingkaran dengan jari-jari r dan titik pusat berada di 
pusat vorteks, yaitu 





Kp V,rao — $ aa dd -y (5.27) 


Ini berarti bahwa kekuatan vorteks, y , sesungguhnya adalah sirkulasi di sekeliling suatu 
lingkaran yang mengurung titik pusat vorteks, dan kelihatannya tidak tergantung pada jari- 
jari lingkaran tersebut. Hasil ini cukup mencengangkan karena aliran vorteks adalah solusi 
persamaan Laplace yang berlaku untuk aliran potensial, dimana seharusnya aliran bersifat 
tak-rotasional, yaitu vortisitas dan sirkulasi seharusnya berharga nol di seluruh medan 
aliran. Sebelum masalah ini kita bahas lebih lanjut secara mendalam, marilah kita periksa 
dulu apakah sirkulasi memang selalu berharga nol untuk aliran elementer lainnya seperti 
angin seragam dan aliran sumber. 

Untuk memudahkan hitungan, sirkit C dimana sirkulasi ingin dihitung dipilih 
berbentuk segi panjang dengan sisi-sisi sejajar sumbu x dan sumbu y, masing-masing 
berukuruan 6x dan 6y. Titik-titik sudut segi panjang tersebut disebut A, B, C dan D 
dimulai dari titik A di pojok kiri bawah dan berputar kebalikan arah putaran jarum jam. 
Sirkulasi pada sirkit ABCD ini adalah 


T - P(udx dy) 


T —AB.u#BC.v -CD.u--DA.v 
T — &x.cos(a) #Sy.sin(a) — &x.cos(a) — 6y.sin(x) — 0 (5.28) 


Jadi aliran angin seragam memang bersifat tak-rotasional. Sirkulasi pada sirkit yang sama 
untuk aliran sumber dapat dihitung sebagai berikut 


B Cc D A 
tr — Judx4 |vay4 Judx 4 | vax 
A B c D 


dimana 
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tj xdx Oo 
udx — — -——In(x 4y 
| 2m! X4y” 41x - 





kel ydy o NA" 
vdy — — -—In(x - 
| Pn PA y An (E 


jadi 








Pn (an) ane) “Ga 


-0 (5.29) 
dengan pengertian bahwa 

RA — VXA? 4 YA? dan seterusnya 
Seperti dapat dilihat aliran sumber memang bersifat tak-rotasional. 

Pembahasan di atas menunjukkan bahwa aliran angin seragam dan aliran sumber 
memang benar-benar bersifat tak-rotasional, sedangkan aliran vorteks agak meragukan. Di 
bawah ini kita akan menunjukkan bahwa aliran vorteks juga benar-benar tak-rotasional 
walaupun persamaan (5.27) seolah-olah tidak membenarkan pernyataan tersebut. Dalam 
pembahasan tentang tornado di bab 3, dijelaskan bahwa dalam mata atau kawasan inti 
tornado, yaitu untuk 0 X r X r, fluida bergerak putar seolah-olah sebagai benda kokoh. Di 
kawasan inti tersebut kecepatan tangensial fluida diberikan oleh rumus berikut 

V-0.r untuk 0srsr, (5.30) 
dimana @ adalah kecepatan putar fluida yang berperilaku seperti benda kokoh. Sirkulasi di 
sekeliling lingkaran dengan jari-jari r £ r, adalah 


T -21.0.r untuk 0srsr, (5.31) 
dan khususnya bila r — r,, maka 
IT, 21.0. (5.32) 


Vortisitas adalah sirkulasi per satuan luas, jadi 


G- aa -20 untuk 0srsr, (5.33) 


Distribusi kecepatan tangensial di luar kawasan inti diberikan oleh rumus berikut 
V—0.1r/r untuk r2»r, (5.34) 
Sirkulasi di sekeliling lingkaran dengan jari-jari lebih besar dari jari-jari inti atau r, adalah 


T -21mr.V, —2m0or” —T, untuk r»r, (5.35) 
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Hasil di atas menunjukkan bahwa sirkulasi pada lingkaran dengan jari-jari lebih besar 
daripada r, itu berharga tetap, sedangkan di dalam kawasan inti, atau untuk r « r,, harga 
sirkulasi adalah fungsi dari jari-jari r. Tetapi sirkulasi di sekeliling vorteks yang diberikan 
oleh persamaan (5.27) itu berharga tetap dan ini berarti bahwa sirkulasi tersebut dihitung 
pada jari-jari yang lebih besar daripada r,. Harga sirkulasi per satuan luas untuk kasus ini 
adalah 


Ni ME (2) (5.36) 


Pi r 





seperti dapat dilihat, persamaan (5.36) menunjukkan bahwa harga C itu tergantung pada 
jari-jari lingkaran r. Ini jelas bertentangan dengan konsep vortisitas di suatu titik yang 
seharusnya bernilai tertentu atau suatu konstanta. Memang vortisitas adalah limit harga 
sirkulasi dibagi luas lingkaran apabila jari-jari lingkaran diambil limitnya mendekati nol. 
Tetapi perlu diingat bahwa dalam konsep aliran vorteks, jari-jari inti yaitu r, telah 
didefinisikan sebagai bernilai nol. Persamaan (5.36) menunjukkan bahwa bila r» 0, 
walaupun r mendekati nol, ternyata harga € adalah nol karena r, — 0 dan @ adalah suatu 
konstanta dengan besaran berhingga. Di titik pusat vorteks sesungguhnya vortisitas tidak 
terdefinsikan karena r, — 0 dan r juga bernilai nol. 

Pembahasan di atas menunjukkan bahwa titik pusat vorteks adalah titik singguler dan 
medan aliran vorteks sesungguhnya bersifat tanpa pusaran atau tak-rotasional kecuali di 
titik singguler itu sendiri. Sekarang marilah kita periksa harga sirkulasi di sirkit sembarang 
yang tidak mengurung titik singguler atau pusat vorteks. Ini harus ditunjukkan selalu 
berharga nol untuk membuktikan bahwa medan aliran vorteks adalah medan aliran 
potensial kecuali di titik singgulernya. Bayangkan suatu sirkit yang terdiri dari 2 lingkaran 
konsentris, dengan jari-jari r, untuk lingkaran luar yang besar dan jari-jari r, untuk 
lingkaran dalam yang lebih kecil. Kedua lingkaran tersebut berpusat di titik pusat vorteks 
dan dua sisi lain dari sirkit itu adalah 2 garis radial pada sudut @, untuk garis radial atas 
dan pada sudut (27—@,) untuk garis radial bawah (lihat gambar (5.7)). Sirkit tersebut 
disebut sirkit ABCD, dimana titik-titik A dan D berada pada lingkaran luar sedangkan B 
dan C berada pada lingkaran dalam. Titik-titik A dan B berada pada garis radial 0,, 
sedangkan C dan D berada pada garis radial 0,,. Sirkulasi di sekeliling sirkit ABCD dapat 
dihitung sebagai berikut 

T - BC. Vo, -CD.0- DA.V,,, #AB.0 


T —(217-0, —0,)r.(y /2m,)—(21-6, — 0, ).r.(y 12115) 
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yaitu 

IT -0 
Jadi ternyata sirkulasi di sekeliling sirkit yang tidak mengurung pusat vorteks ternyata 
berharga nol. Dapat diamati bahwa hasil tersebut berlaku umum untuk semua harga 0,, 
O,, r, dan r,, dengan catatan bahwa 0 « (0, t0,) «21 dan 1,22. Sirkit ABCD 
dengan harga 9, — 0, — 0, r, —r, dan r » r, dapat dilihat pada gambar (5.7b). Bila jari- 
jari inti vorteks, r,, sekarang diambil limitnya mendekati nol, maka akhirnya sirkit ABCD 
menjadi seperti ditunjukkan pada gambar (5.7c). Sirkit terakhir ini sesungguhnya adalah 
seluruh kawasan di dalam lingkaran denga jari-jari r,, berharga sembarang, tetapi tidak 
mencakup titik pusat vorteks yang merupakan titik singguler. 

Pembahasan di atas menunjukkan bahwa medan aliran vorteks adalah medan aliran 
potensial atau tak-rotasional, kecuali di titik singgulernya. Perhatikan bahwa dalam 
gambar (5.7c), lingkaran yang mengurung titik singguler telah dipotong sepanjang suatu 
garis radial. Walaupun garis potong itu kelihatannya hanya ada satu, tetapi seandainya 
lingkaran itu dipotong dengan gunting maka sesungguhnya ada dua garis potong, yaitu 
garis potong atas, AB, dan garis potong bawah, CD. Dipandang dari segi teori variabel 
kompleks, keberadaan garis potong tersebut adalah untuk memenuhi persyaratan adanya 
turunan dari fungsi kompleks, f' (z). Persyaratan tersebut adalah bahwa fungsi f(z) harus 
bernilai tunggal (single-valued). Kita tahu bahwa 1 titik pada bidang kompleks z itu dapat 
mewakili sejumlah tak berhingga bilangan kompleks. Misalnya saja semua bilangan 
kompleks yang diwakili oleh rumus 

z— r.exp|i(0 t2nn)| Mm-0,1,2,... 
dalam kenyataannya diwakili oleh 1 titik yang sama. Dengan adanya garis potong AB 
(atau CD) kita telah memaksa n hanya boleh berharga n — 0 saja, dan dengan demikian 
kita memaksa fungsi f(z) bernilai tunggal. Kenyataan bahwa aliran vorteks adalah aliran 
potensial diseluruh medan aliran kecuali dititik pusatnya itu sangat penting dipandang dari 
segi aerodinamika pesawat terbang. Teori Kutta-Joukowski menyatakan bahwa gaya 
angkat yang beraksi pada suatu benda itu berbanding lurus dengan sirkulasi disekeliling 
sirkit yang mengurung benda tersebut. Teori aliran potensial itu didasarkan pada anggapan 
bahwa aliran udara bersifat tak-rotasional, yaitu sirkulasi di sekelilling sirkit apapun 
dimedan aliran itu harus berharga nol. Bila sirkulasi selalu berharga nol maka teori aliran 
potensial tidak akan bermanfaat untuk menghitung gaya angkat. Kenyataan bahwa titik 
pusat vorteks bersifat rotasional, dan sirkulasi disekeliling sirkit yang mengurung titik 
singguler tersebut tidak bernilai nol, berarti bahwa teori aliran potensial dapat 
dimanfaatkan untuk menghitung gaya angkat. 


59 


Suatu hal lain yang perlu diamati adalah bahwa titik vorteks itu sesungguhnya adalah 
titik potong suatu garis vorteks dengan panjang tak-berhingga, sejajar sumbu-z, dengan 
bidang datar x-y. Berbeda dari aliran sumber yang dapat bersifat 2-dimensional sebagai 
garis sumber, atau 3-dimensional sebagai titik sumber, aliran vorteks ternyata hanya dapat 
diwakili oleh garis vorteks dan tidak dapat berbentuk titik vorteks. Tetapi perlu disadari 
bahwa garis vorteks itu boleh berbentuk kurva lengkung atau garis lurus. Ini berbeda dari 
garis sumber yang hanya dapat berbentuk garis lurus. 

Aliran elementer keempat yang akan dibahas selanjutnya adalah aliran doblet yang 
fungsi potensial kompleksnya diberikan oleh rumus (5.19). Fungsi potensial kecepatan dan 
fungsi arus untuk kasus ini adalah 











p x 
9 y) 5 —.G—T 
21 xXx 4y (537) 
PENA AN 
w(x,y) - In Pa 
Komponen kecepatan untuk aliran doblet adalah 
Sesaln AA 
up 21 (4 
yu 2xy 
.--——— 5.38 
v(x,y) an May (5.38) 
Komponen kecepatan dalam koordinat polar adalah 
H cos@ 
V, (1,90) ———. 
(7:8) 217 1 
NA an PAR Na (5.39) 
At 21 r 3 


Rumus matematis untuk garis arus adalah 
ey )ay-0 
In 
yang dapat ditulis ulang menjadi 


2 2 
2 afve h | ) (5.40) 
Ay, 41, 


Persamaan diatas menunjukkan bahwa garis arus aliran doblet adalah lingkaran-lingkaran 





dengan jari-jari rp dan titik pusatnya berada dititik (O,trk), dimana 
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ni — (5-41) 





7. AT, 
Rumus matematis untuk kurva isopotensial adalah 
2 2 
Ke — apa — (5.42) 
AT, 4, 


Jadi kurva isopotensial berbentuk lingkaran berjari-jari rk dan titik pusatnya berada di (Ht 
Tk:0). Untuk yu berharga positif arah aliran pada lingkaran diatas sumbu-x adalah searah 
putaran jarum jam, sedangkan pada lingkaran dibawah sumbu-x adalah kebalikan arah 
putaran jarum jam (lihat gambar (5.8)). 

Sekarang bayangkan suatu lingkaran berjari-jari r dan titik pusatnya berada dititik 
doblet. Fluks massa dan sirkulasi pada lingkaran tersebut dapat dihitung sebagai berikut 


v-fv.rao-4( 5. Jeostopao -0 (5.43) 





r-hu.rao ff 3.5 Jsintopad -0 (5-44) 


Seperti dapat dilihat fluks massa pada lingkaran r ternyata berharga nol, mirip dengan 
kondisi untuk aliran vorteks. Tetapi aliran doblet berbeda dari aliran vorteks karena 
sirkulasi pada lingkaran r ternyata berhargan nol untuk kasus aliran doblet. 

Pengertian yang lebih mendalam tentang aliran doblet dapat dipelajari sebagai 
berikut. 

Bayangkan suatu garis sumber dengan kekuatan (0/21) diletakkan menembus 
bidang x-y pada titik (ta,0). Kemudian suatu garis serap dengan kekuatan (6/27) 
diletakkan pada titik (-a,0) dan kita pelajari pola aliran yang dihasilkan. Garis sumber yang 
diletakkan di (ta,0) akan menimbulkan suatu medan aliran dengan potensial kompleks 
berikut 





D,(z—a)--“Infz—a) (5.45) 
21 
Sedangkan potensial dari garis serap di (-a,0) adalah 
D,(zta)-— Infzta) (5.46) 
21x 
Potensial kompleks dari kombinasi linier sumber dan serap tersebut adalah 
0-0, 10, -- | 272) (5.47) 
21 zta 
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Bila dirumuskan bahwa 

2 5Z-—a-—r.exp(iO)—a—r,.exp(iO,) 

z, 5zta—r.exp(iO) ta —r,.exp(iO,) (5.48) 
dimana 

r 5 ta? —2a.r.cos9 

rr ta? #2a.r.cos0 

tan0, —r.sin@/ (rcosO—a) 

tan0, —r.sin0/ (rcosO--a) (5.49) 
maka fungsi potensial (5.47) dapat diuraikan menjadi 


2 


D(z) - 2 (5) cxp(i(0, -)| (5.50) 


atau 


6(T,0) Hiw(r,0) — Sin L yi (0, —0,) (5.51) 
21 r, 2x 


Dengan sedikit pengolahan dapat ditunjukan bahwa 





00 | Kay) 
b(x,y) Lah La (5.52) 
Va) tan (5.53) 


Kurva isopotensial diberikan oleh rumus 
(x—a)P ty' — ka tx ta y) 
k, - exp(4nd, /0) (5.54) 
Persamaan (5.54) dapat diolah lebih lanjut menjadi 
(Ary na? 
r, Sa(k, #1)/(k, —1) (5.55) 
Apabila harga a sekarang diambil limitnya mendekati nol, tetapi k,, dipilih sedemikian 
rupa sehingga rp berharga tetap, maka persamaan (5.55) dapat diolah menjadi 
(ant 5 (5.56) 
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Perhatikan bahwa persamaan (5.56) itu tepat sama dengan persamaan (5.42) yaitu 
persamaan untuk kurva isopotensial aliran doblet. Kurva garis arus diberikan oleh rumus 
(5.53) yang dapat ditulis menjadi bentuk berikut 
Ryan) 5 ta? 
r--a/k, 
k, -tan(27y, /o) (5.57) 

Perbandingan dengan persamaan (5.40) menunjukkan bahwa persamaan (5.57) adalah 
rumus untuk garis arus aliran doblet. Hasil diatas menunjukkan bahwa aliran doblet 
ternyata dapat diperoleh sebagai hasil kombinasi linier sebuah garis sumber dengan sebuah 
garis serap, yang diletakkan di titik yang sama. Walaupun kekuatan sumber dan serap itu 
sama (tetapi berlawanan tanda), untuk memperoleh aliran doblet harus dibuat anggapan 
bahwa sumber dan serap tersebut tidak saling meniadakan. Perkataan doblet berarti 
sepasang singgularitas dan kita telah membuktikan bahwa doblet sesungguhnya adalah 
pasangan sumber dan serap dalam limit tertentu. Dibawah ini akan dibuktikan bahwa 
doblet dapat juga diartikan sebagai pasangan vorteks positif dan negatif di titik yang sama, 
dalam limit tertentu. Fungsi potensial kompleks dari kombinasi linier vorteks positif di titik 
(O,ta) dan vorteks negatif dititik (O,-a), masing-masing dengan kekuatan (y/21), adalah 

Oz) ——i-Infz-ia) tiD-Infzia) 

21 21 

Persamaan diatas dapat ditulis ulang sebagai 

9(z)- Lal 2) (5.58) 

21x z—ia 

Bila dirumuskan bahwa 

px 4(ytay 

rx #(y-a) 

tan(O,) — (y#-a)/x 

tan(0,) — (y-—a)/x 

tan(O, —0,) - 2ax/(x'4y'—2”) (5.59) 
maka persamaan (5.58) dapat diuraikan menjadi 





&(z) - - (0, -0,)-i 1 Te 


Ty 


atau 
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desy) — tan 2 ) (5.60) 


2x Kya? 
2 2 
Y x #(yta) 
»y) 5 —In —— .61 
y(x,y) am | 020) Be) 
Garis-garis arus diberikan oleh persamaan berikut 
(yen) —r —a, (5.62) 
dimana 
Ksa(1t-k,)/(1-—k 
k5a(1tk3)/(1—k3) (5.63) 


ks —exp(4ny, /y) 
Bila a diambil limitnya menjadi nol tetapi rk berharga tetap maka persamaan (5.62) 
disederhanakan menjadi 

(yr) (5.64) 
yang juga merupakan persamaan untuk garis arus aliran doblet (lihat persamaan (5.40)). 
Dengan cara yang sama dapat juga ditunjukkan bahwa persamaan (5.60) dapat diolah 
menjadi persamaan atur untuk kurva isopotensial aliran doblet. 

Persamaan Laplace adalah persamaan diferensial derajat-2 yang bersifat linier. Ini 
berarti bahwa suatu solusi yang lebih rumit selalu dapat diperoleh sebagai hasil kombinasi 
linier 2 atau lebih solusi-solusi yang bersifat lebih elementer. Sebagai contoh kita telah 
menunjukkan bahwa kombinasi linier aliran sumber dan serap, atau vorteks positif dan 
negatif ternyata dapat menghasilkan solusi baru, yaitu doblet. Di bawah ini kita akan 
mempelajari pola aliran yang dihasilkan apabila aliran angin seragam digabungkan dengan 
aliran doblet. Fungsi potensial kompleks untuk aliran ini adalah 


02) 11 (5.65) 
21 z 
atau 
2 
D(z)- st 
z 
Re—- (5.66) 
2x 


Mengingat bahwa z.z — le, maka persamaan diatas dapat diolah lebih lanjut menjadi 
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